DISTRIBUTION OF THE FIRST PARTICLE IN DISCRETE ORTHOGONAL 

POLYNOMIAL ENSEMBLES 



ALEXEI BORODIN AND DMITRIY BOYARCHENKO 



Abstract. We show that the distribution function of the first particle in a discrete orthogonal polynomial 
ensemble can be obtained through a certain recurrence procedure, if the (difference or q-) log-derivative of the 
weight function is rational. In a number of classical special cases the recurrence procedure is equivalent to the 
difference and q-Painleve equations of ]lo[ , [^?|. 

Our approach is based on the formalism of discrete integrable operators and discrete Riemann-Hilbert problems 
developed in Hj, jij. 
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1. Introduction 

1.1. The basic problem considered in this paper is the following. Let X be a locally finite subset of M. and 
w : X — > M>o be a positive-valued function on X with finite moments: 

\x\ n w(x) < oo, n = 0, 1, . . . . 

Fix a positive integer k (the number of particles) and consider the probability measure on all fc-point subsets 
of X given by 

k 

Prob{xi, ... ,Xk} = const • 1Q (x % - x ) 2 w(xi). (1.1) 

\<i<j<k i=l 

We are interested in the distribution of max{ii, . . . , Xk} with respect to this measure. 

The problem is motivated by random matrix theory on one side, and by combinatorial and representation 
theoretic models on the other one. 
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In random matrix theory, probability measures of the form 

k 

const • Y\ ( x i ~ x j) 2 W w(xi)dxi 

l<i<j<k i=l 

on fc-point subsets of R, with w{x) being a smooth function on a subinterval of R, play a prominent role. Most 
computations for such models are conveniently done by means of the orthogonal polynomials associated with 
the weight function w(x). On this ground, these measures are often called orthogonal polynomial ensembles. See 
Jl4) , [lq] and references therein for a further discussion. 

The problem of describing the distribution of the max{xi} in the continuous setting for the classical weights 
has been solved in the following sense: the distribution function was explicitly written in terms of a specific 
solution of one of the six (2nd order nonlinear ordinary differential) Painleve equations. It was done in pOfl for 
the Hermite weight w(x) — exp(— x 2 ), id, and for the Laguerre weight x a exp(-x), x > 0; in pQ ), H for the 
Jacobi weight (1 - x) a (l + x) b , x e (— 1, 1); and in [Q, || for the quasi- Jacobi weight (1 - ix) s (l + ix) s , x G R.Q 
Thus, it is natural to ask what would be an analog of these results when we take w to be a classical discrete 
weight function. 



On the other hand, in recent years the random variable maxjxi} with Xj's distributed according to (1.1) 
with certain specific weights, came up as the main quantity of interest in a number of problems originating in 
combinatorics, first-passage percolation, representation theory, and growth processes, see e.g. 0, Q, [jij), || 
and references therein. 

1.2. In order to state our first result we need to introduce more notation. Let us denote the points of X by 7r s , 
s = 0, 1, . . . , N, with 7r < tti < ■ ■ ■ or 7r > 7Ti > • • • . Here N = |3£| — 1 may be finite or infinite. We use the 
following two basic assumptions: 

• There exists an affine transformation a : R — > R such that aTi s+ i — ir s for all s, < s < N. 

• There exist polynomials P(x) and Q(x) such that 

w( - w ;-^ = ^4, i < s < n, 

w(tt s ) Q{n s ) 

and P(n ) = 0. 

The orthogonality data for a number (but not all) hypergeometric polynomials of the Askey scheme satisfy these 
assumptions, see §^ below for details. f] 

We prove that under the two conditions above, there exists a certain recurrence procedure which computes 
the gap probability 

I Prob{max{:Ei} < 7r s }, if ttq < 7Ti < • • • , 



D s = Prob {xi <£ {7r s , 7r s+ i, . . . } for all i} 



|Prob{min{a;i} > ir s }, if ttq > m > ■ ■ ■ 



with x^s distributed according to (1.1). In fact, the recurrence procedure produces a sequence (A s , M S (C)), where 



A s is a nilpotent 2 by 2 matrix and M S (Q is a matrix polynomial 

M S (C) = M { pC l + ■■■ + M s (0) , e Mat(2, C), 

of degree / = maxjdegP, degQ}. The elementary step of the recurrence is provided by the equality 

i + -r*—) M *(0 = M a+1 (o (i + y^ 2 -) ■ (i.2) 



1 Continuous problems of this type have been extensively studied. We refer to the introduction of for a brief review and 
references. 

2 For some classical families of polynomials both assumptions are satisfied but the orthogonality set X is not locally finite. We 
extend our results to those cases, see 9d below. 
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It is not hard to see that if det M s (tt s +i) ^ (which is always the case in our setting) then (1.5) defines 
(A s+ i, M s +i) uniquely provided that we know {A s ,M a ). However, the existence of {A s +\, M s +i) is not obvious 
and needs to be proved.^ Again, in our setting it always holds. 
We then show that the ratio 

-l 
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D s+ 2 
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s+l 



and (A s+1 ,M t 



(0 

S+l! ' 



Ml) 



is an explicit rational function of (A s , Ms 

Since D s — Prob{max{xi} < 7r s } is nonzero only if s > k (recall that k is the number of x^s in (1.1)), it 
is enough to provide the initial conditions -D&, -Dfc+i, -Dfc+2, A^, M/.(£) in order to be able to compute D s for 
arbitrary s. These initial conditions are readily expressed in ter ms of {n s } and {w(tt s )}, see §[| below. 

For certain classical weights w the recurrence relation (1.2) can be substantially simplified. To illustrate 
the situation, let us consider X = Z>o and w(x) = a x /x\, where a > is a parameter. This weight function 
corresponds to the Charlier orthogonal polynomials. 

In this case, A s and M s (() can be parameterized by three scalar sequences a s , b s , c s as follows: 



M s (0 



c + 



b s 
a/c s 



b s c s 
a 



As 



(k 



-1 
l/{a s c s 



-a s c s 
1 



Then the equality (1.2) leads to the following recurrence relations: 

(ps + aa s )(fc + bs + aa s ) 
aa s (s + 1 + b s + aa s ) 



1 



5s+l 



1 - a s +i 

Cs+l 



(s + 1 + k + b s + aa s ), 



k + b s + aa s 

The connection of these sequences and the distribution D s is given by 



(1.3) 

(1.4) 

(1.5) 
(1.6) 
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s+l 
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l s+l) 



a(s + 2)a 2 s+1 c s +i 



The corresponding initial conditions can be found in subsection 
Under the change of variables 

fa = a^ 1 , 



11.2 



(1.4)-(1.5) turn into 



fsf: 



s+l 



ag s 



9s + 9s+i 



(g s — s— l)(g s + k - s - 1)' 
a s + l 



k + 2s + 3. 



/s+i 1 - /s+l 

This recurrence is immediately identified with the difference Painleve IV equation (dPIV) of |l7| . 

1.3. It turns o ut that the situation for the Charlier weight described above is rather typical. We are also 
able to reduce (L2) to scalar rational recurrence relations for the weight functions corresponding to Meixner, 
Krawtchouk, q-Charlier, alternative q-Charlier, little q-Laguerre/Wall, little q-Jacobi, and q-Krawtchouk orthog- 
onal polynomials. In the appropriate variables, Meixner and Krawtchouk cases lead to dPV of [[l7|, little q-Jacobi 
and q-Krawtchouk lead to q-P VI of jl0| , |jl7) , and q-Charlier and little q-Laguerre yield a certain degeneration 
of q-PVI. 



3 In fact, for 2 by 2 matrices one can easily see that (j4 s _|_i, M 



(0 

s + l' ' 



, M^9j) are rational functions of (A B , 



M 0) ). 
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It is remarkable that in almost all the cases we can solve explicitly, we end up with one of the equations of 
Sakai's hierarchy which was constructed out of purely algebraic geometric considerations, see |l7| . (We were not 
able to see such a reduction in the alternative q-Charlier case, but we do not claim that there is none.) So far 
we have not found a conceptual explanation for this fact. 



One can notice, however, that recurrence relations originating from (1.2) must have some kind of singularity 
confinement property. (This property was the starting point of Sakai's work.) For example, the parameterization 
( |l.3| ) does not make much sense if, say, A s has a zero (2,1) element. Then the values of a s and c s are not 
well-defined. In terms of the recurrence relations, this is reflected by vanishing of one of the denominators 
in ( |1.4| )-( |L6| ). However, the matrix sequence {(A s , M s )} does not feel this singularity, which means that the 
sequences {a s }, {b s }, {c s } can be "continued through" their singular values. Of course, all Sakai's equations have 
this kind of singularity confinement by construction. 

Let us also point out that it is not clear at this point whether the weights of higher hypergeometric polynomials 
of the Askey scheme will also le ad t o one of Sakai's equations. All the cases that we were able to solve explicitly 



have linear matrices M s (£) in (1.2), while, say, for the Hahn weight M S (C) is quadratic. Handling such cases 
seems to be a problem of the next level of difficulty. It remains an interesting open problem to derive explicit 
rational recurrence relations for degM, = 2. 

For Charlier and Meixner weights, D s can also be written as Toeplitz determinants with symbols 

(1 + z) fe exp(az- 1 ) and (1 + z) k (l + bz^f (1.7) 

respectively, see [pH , ||. Here a, b, c are parameters. Among previous results on the subject let us mention 

• the derivation of dPII for Toeplitz determinants with the symbol exp(8(z + z -1 )), see 0, [Q, jlj (note also 
the derivation of the same equation for the quantity closely related to these Toeplitz determinants in p6| , 

©); 

• derivation of dPV for Toeplitz determinants with symbol (1 + z) k (l + bz 1 ) c and k being not necessarily 
integral in Q ; 

• derivation of rational recurrence relations for Toeplitz determinants with symbols of the form 

exp(P 1 (z) + P 2 (z- 1 )) _ d lZ y'i(l - d 2 z)<{l - d^z-Y 1 '^ - d^z- 1 ) 1 * (1.8) 
where P\ and P 2 are polynomials with jdegPi — degP2| < 1, 7i ? 72 i 7i' ; 7a i dt, d% are constants, see [jjj 



Interestingly enough, for the symbols ( |l.7|) , which are special cases of (1.8), the relations of [y do not seem 
to have much in common with those of[Q and the present paper. 

1.4. The methods used in this paper are based on the formalism of discrete integrable operators and discrete 
Riemann-Hilbert problem (DRHP) developed in ||, j|. The first step is to represent D s as a Fredholm deter- 
minant of an integrable operator: D s — det(l — K s ), where K s is an operator in £ 2 ({ir s , 7r s+ i, . . . }) with the 
kernel 

*.(*, y) = j, — p — p^)p^y)-p^)p^y) y^%y . 

\\Pk-i\\ £H3£tW) x-y 

Here pk and pk-i are monic kth and (fc — l)st orthogonal polynomials on X with respect to the weight function 
w. Using the results of the computation of such a Fredholm determinant can be reduced to solving a 

DRHP on {ttq, . . . , 7r s _i} with a jump matrix easily expressible in terms of w. Our assumptions on X and w, see 
above, then allow us to obtain a Lax pair for the solution m s (C) of this DRHP, which has the form 

m s+ i(C) = (l+ m ^), m °( a = MsiOms+iiOD- 1 ^)- (1-9) 

Here M S (Q is a matrix polynomial and -D(C) is a fixed diagonal matrix polynomial. Compatibility condition for 



this pair of equations is exactly (1.2). 

The paper is organized as follows. In we reduce the problem of computing discrete orthogonal polynomials 



with a given weight to a DRHP. In §|| we derive the Lax pair (l.S). In §^ we show how to express the Fredholm 
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determinants D s in terms of (A s , M s ). In §|] we prove that the compatibility condition ( |l.2| ) always has a unique 
solution. In §[j]we derive the initial conditions Ak, Mk, P>k, -Dfc+i- In ^ we write down explicitly the Lax pairs for 
14 families of discrete hypergeometric orthogonal polynomials of the Askey scheme. In §||we solve (1.2) in terms 
of rational recurrences for the matrix elements of A s and M s for deg M s = 1. In §p|we show how to reduce (1.2) 
to dPIV and dPV equations in the case when degM s = 1 and a( — ( — 1. In § |To| we reduce (1.2) to q-PVI or its 
degeneration in the case when degM, = 1 and a( — g ±1 C- Finally, in §[TT| we solve (1.2) in terms of difference and 
q-Painlcve equations for 7 families of classical orthogonal polynomials. At the end of the paper we also provide 
a few plots of the "density function" (difference or g-derivative of D s ) obtained using the formulas of § p"l| . 

This research was partially conducted during the period the first author (A. B.) served as a Clay Mathematics 
Institute Long-Term Prize Fellow. He was also partially supported by the NSF grant DMS-9729992. 

1.5. The following notation is used throughout our paper. For an integer k, we write Z>^ = {k, k + 1, fc+2, . . . }. 
If a,q £ C and k £ Z>o, one defines the Pochhammer symbol and its g-analogue (often also called the q-shifted 
factorial) by 

(a) ■= 1, (a)* := a(a + l)(a + 2) • • • (a + k - 1) if k > 1 

and 

(a; g) := 1, (a; q) k := (1 - o)(l - aq)(l - aq 2 ) ■■■{!- aq^ 1 ) if k > 1, 
respectively. One usually writes 

r r 

(ai,...,a r )k = Y\_(aj)k and (a%, . . . , a r ; q) k = JJ(aj-; q)k- 

.7=1 J'=l 

If r, s £ Z>o and a\, . . . , a ri b\, . . . , 6 S , z, g £ C, the hypergeometric series and the 6asic hypergeometric series are 
defined by 

(ai, . . .,a r )k z k 



p 



oi, . . . , a r 
bi,...,b s 



E 

fe=0 



(&i, . . .,b s ) k k\ 



and 



respectively. 



7<t>i 



( at,. 


. . ,a r 






..,b s 





(,ai,...,aT-;gJfc ^ 1 x(i +s - r )fc g (i+ a -r)(S) 2 
^ (bi,...,b s ;q) k (q;q)k' 



k=0 



2. Discrete Riemann-Hilbert Problems and Orthogonal Polynomials 

2.1. In this section we explain how solutions of discrete Riemann-Hilbert problems (DRHP) for jump matrices 
of a special type can be expressed in terms of the corresponding orthogonal polynomials. Let X be a discrete 
locally finite subset of C, and let w : X — » Mat(2, C) be a function. As in H, j|, we say that an analytic function 

m : C\X — ► Mat(2,C) 

solves the DRHP (X, w) if m has simple poles at the points of X and its residues at these points are given by the 
jump (or residue) condition 

Resm(C) = lim(m(C)w(x)), x £ X. (2.1) 

£ — X £ — >x 

Lemma 2.1. If m(£) is a solution of the DRHP (X, w) and the matrix w(x) is nilpotent for all x £ X, then the 
function detm(C) is entire. If, in addition, detm(C) — > 1 as ( -> oo, then detm(C) = 1- 

Proof. For each x £ X, the jump condition (2.1) implies that the function m(Q ■ (i — (C — x)~ 1 w(x)) is analytic 
in a neighborhood of x. Since w{x) is nilpotent, this product has the same determinant as m(C), which shows 
that det77i(C) has no pole at x. The second statement of the lemma follows from Liouville's theorem. □ 
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2.2. We now assume that the matrix w(x) has the following form: 

uj(x) 



»(*) = ^ y), (2.2) 

where u> : X — > C is a function. Recall that a collection {P n (C)}5£L Q of complex polynomials is called the collection 
of orthogonal polynomials associated to the weight function uj if 

• P n is a polynomial of degree n for all n = 1, 2, . . . , and Po = const; 

• if ?7i ^ n, then 

^P m (a;)P„(zV(a;) =0. 
We will always take P n to be monic: P n (x) — x n + lower terms . 

In order for the definition to make sense, we assume that all moments of the weight function to are finite, i.e., 

the series \u!(x)x j \ converges for all j > 0. (2.3) 

Under this condition, one can consider the following inner product on the space C[C] of all complex polynomials: 

(/(C), 5(0)- := E ffr)9(*M*)- 

It is clear that there exists a collection of orthogonal polynomials {P n (C)} associated to uj such that (P„, P n )u ^ 
for all n if and only if the restriction of (•, -) w to the space C[£]- d of polynomials of degree at most d is 
nondegenerate for all d > 0. If this condition holds, we say that the weight function uj is nondegenerate. In this 
case, it is also clear that the collection {P n } is unique. 

Remark 2.2. If the set X is finite, one has to modify the definitions above. Indeed, if X consists of N + 1 points 
(N G Z>o), the inner product (•, -) u is necessarily degenerate on C[£]- d for all d > N. So, instead, we require 
that (-, -) u be nondegenerate on C[£]- d for < d < N, and we are only interested in a collection {P n (0}n=o 01 



orthogonal polynomials of degrees up to N. On the other hand, the condition (2.3) is empty in this case 



Remark 2.3. If the values of the weight function uj are real and strictly positive and the orthogonality set X is 
contained in R, then u> is automatically nondegenerate. 

2.3. The first basic result of the paper is 

Theorem 2.4 (Solution of DRHP). Let X be a discrete locally finite subset o/C and lu : X — > C be a nondegener- 



ate weight function satisfying (2.S). Let {P n (C)}n=a be the collection of monic orthogonal polynomials associated 



to uj, where N = card(X) — 1 G Z>o U {oo}. Assume that the jump matrix w(x) is given by (2.i ). Then for any 



k = 1, 2, . . . , N, the DRHP (X, w) has a unique solutions mx(C) satisfying the asymptotic condition 



m*(CH C n % ) =I + o(-A asC-^, (2.4) 



where L is the identity matrix. Lf we write 



then m£(Q = P k (0 and mf{() = (P k ^, Pfe-r)" 1 • P fe -i(C). 

Remark 2.5. The asymptotic condition (|2.4|) needs to be made more precise. Indeed, if X is infinite, then the 



LHS of (2.4) has poles accumulating at infinity. In this case, we require that the asymptotics be uniform on a 
sequence of expanding contours (e.g., circles whose radii tend to oo) whose distance from X remains bounded 
away from zero. A similar remark applies to all asymptotic formulas below. 
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Remark 2.6. A result for continuous weight functions similar to Theorem [2.4| was proved in 
Proof. Fix a natural number k < N and define a matrix-valued function 

^ 12 \ 

22 : C — > Mat(2,C) 



™ 12 (C) = E 



by 

™ n (0 = P k (0, 

P k (x)u>(x) 
C-x ' 

xGX * 

m 21 (C) = cPn((), 

22,, n \ - gfe-l(gMg) 

161 S 

where c is the unique constant for which m 22 (0 = C~ + 0(C _fc_1 ) as £ — > oo (we will show below that such a c 
exists, and in fact c = (Pfc-i, Pfc-i)^ 1 ). We first observe that ( [2. 3D implies that the series for m 12 (0 converges 
uniformly and absolutely on compact subsets of C \ X, and hence m 12 (£) is analytic on the complement of X. 
Moreover, since for each id, the series 



E 



P k (y)uj(y) 



C-y 

v ex\{x] s y 

converges uniformly and absolutely on compact subsets of (C \ X) U {x}, we see that m 12 (0 has a simple pole at 
x, with the residue given by 

Resm 12 (C) = P k (x)ui(x). 

Similarly, m 22 (0 is analytic away from X, with a simple pole at each i£l satisfying 

Resm 22 (C) =c-P k -i(x)w(x). 



This shows that our matrix m(0 satisfies the jump condition (2.1). To verify the asymptotic condition (and find 
the constant c), note that 

m n (0-r fc = 1 + 0(1/0, m ai (C)-C* = 0(VC). asC-oo. 

Next, we write 

By the definition of orthogonal polynomials, we have 

E Pk{x)uj{x)x l = 0, < i < Jfe — 1, 

x£X 

and hence substituting ( |2.5| ) into the definition of m 12 (0 yields 

m 12 (0 = OiC^ 1 ) as C -> oo. 



Also, fl2.5|) gives 

m 22 (0 = c • E • C k + 0(C fc_1 ) as C -» oo. 
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Hence, if we set 

\x£X J 

then the matrix m(Q satisfies all the conditions of the theorem. 

To prove uniqueness, let mx be an arbitrary solution of the DRHP (X, w) satisfying the asymptotic condition 
( |2.4| ). Since the functions m and mx have the same (simple) poles and satisfy the same residue conditions at 
these poles, it is clear that the function /(£) = m(Q~ 1 • mx(Q is entire (note that m(£) is invertible by Lemma 
|2.ip . The asymptotic conditions on the matrices m and mx imply that f(() — » / as £ — > oo. By Liouville's 
theorem, f = I, which means that m = mx- □ 

3. Lax pairs for solutions of DRHP 

3.1. Let X be a discrete locally finite subset of C, let w : X — ► Mat(2,C) be a function, and fix a natural 
number k < card(X). If w arises from a nondegenerate weight function w on I with finite moments, as in §^|, 



then from Theorem 2.4 we know that the DRHP (X,u>) admits a unique solution mj(() with the asymptotics of 
diag(C fc , C~ fc ) at infinity. 

Convention. From now on, we assume that the weight function under consideration is everywhere real and 
strictly positive, and the orthogonality set is contained in R. 

If 3 C X is any subset of cardinality > k, it follows from Remark |2.3| and our convention that the restriction 
of u to 3 is also nondegenerate, whence by Theorem |2.4|, the DRHP (3, w\S) has a unique solution m^(Q such 



that™ 3 (0- ( C ~*<°) ^asC 



Let us now assume that the set X is parameterized as X — {ir x }^_ , where N £ Z>o U {oo}. For every s > 0, 
we consider the subset 3s := {^x}x<s-i Q X. If s > k, then card(3s) > k, so by the previous paragraph, we 
have the corresponding solution m s (() :— m^ B (Q. It can also be shown that even though card(3fc) = k, the 
DRHP (3fc, mL ) still has a unique solution mk{Q- Indeed, uniqueness can be proved by the same argument 



as in Theorem |2.4] , and in Proposition OA below we give an explicit formula for mk(C)- Note also that m s (C) 1S 
defined for all s > k; in particular, if N is finite, we have m s (Q = mx(C) for all s > N. 

Our next basic assumption is: 

there exists an affine transformation a : C — > C such that an x+ i — ir x for all < x < N. (3-1) 

The Lax pair in our setting will consist of two equations, one of which relates m s +i(C) with m s (£) and the other 
one relates m s (aC,) with m s +x(C)- We will denote the derivative of a (which is a constant) by m so that 

<t(Ci) - <KG0 = V ■ (Ci - Ca) forallCi,C 2 eC. (3.2) 

The cases of special interest are those when X is the orthogonality set for one of the families of discrete hyperge- 
ometric orthogonal polynomials of the Askey scheme. In this situation, X is a subset of either a one-dimensional 
lattice or a one-dimensional (/-lattice in C, and a is given by either cr£ = £ — 1 or <r£ = q^ 1 ^. These cases will be 
treated in greater detail in §(?]; for now we concentrate on the general theory. 

3.2. The main result of this section is 

Theorem 3.1 (Lax pair). For each s = k, k + 1, k + 2, . . . , let 3s = {^ x G ■£ \x < s — 1} , and let m s (£) be the 
unique solution of the DRHP (3 S) toL ) smc/i £/ia£ m s (() • ( ^ ° fc J — > I as ( — * oo, where w is given 



|3> , „„„„ ,.„ svsy ^ o ^ 

w(x) 




w(ir x ) 



4 To simplify notation, we assume that the weight function w(x) is always defined for x £ Z>q, x < N. 
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(a) For each s 6 1^>k, s < N, there exists a constant nilpotent matrix A s such that 



m s+1 (0 = (l+ m,(C). (3.3) 

(b) Assume that there exist entire functions gJi(C), d 2 (Q such that 

i^^^.fjlM fcr all 1 < * < JV, 
w[x) d 2 {n x ) 

di(no) — Oj an d d 2 (a~ ttn) = if TV is finite. 

Let 

Then for every s € Z>fc, we /lai/e 

m s K) = M.(C)m. + i(C)D- 1 (0. (3.4) 

where M s (£) is entire^. 

(c) With the assumptions of part (b), suppose that the functions di(£) and d 2 (Q are polynomials of degree at most 
n in £, and wriie di(C) = AiC" + (lower terms), e?2(C) = A2C™ + (lower terms). 5ei ki = ?7 Ai, K2 = vT k \ 2 . 
Then the matrix M s (£) is polynomial of degree at most n in £, with the coefficient of Q n equal to diag(«i , k 2 ) ■ 



Equations (3.2) and (3.4) constitute the Lax pair. 

Remark 3.2. As follows from the proof below, the condition d 2 (a~ 1 TTpf) = in part (b) of the theorem is only 
required to assert that M s (£) is entire for s > N. If d^a^ 1 n m) ^= 0, then (b) and (c) still hold for all s, k < s < N. 

Proof, (a) Fix s G Z>o, s < N, and consider the matrix-valued function N(() := rn s+ i(C)m j r 1 (C) (recall that 
m s (C) is invertible by Lemma [O]). It is clear that N(() has only one simple pole, at £ = 7r s . Hence the function 
-W(C) — (C ~ 7Ts) _1 ^4s is entire, where A s = Resf =7Ts -/V(C). But m s (C) and m s +i(C) have the same asymptotics at 
i nfin ity, so JV(£) — (f — 7r s ) _1 A s — > / as £ — > 00. By Liouville's theo rem , N(Q — (£ — 7r s ) _1 yl s = /, which gives 
( |3.3|) . Taking determinants of both sides of ( |3.3|) and using Lemma |2.l| gives det(/ + (C - tts^As) = 1 for all 
£, which forces the matrix A s to be nilpotent. 

(b) We have M s (£) = m s (<7C)-D(C)m7+i(C)- Therefore it is clear that M S (Q is analytic away from {n , tti, . . . ,n s }. 
Let 1 < x < s. Then for £ near tt^, we can write 

m s (a0=H 1 (0(l+ W ( X ~V ) and mjUO = (i - H 2 (Q, 

where H\ and i?2 are analytic and invertible matrix- valued functions defined in a neighborhood of tt x . So M s (£) 
is analytic near if and only if so is the product 

i.e., if and only if 



u)(x) d 2 (n x ) 
Similarly, we note that m s (cr£) is analytic near ttq, whereas 

-«<o = (^)«<o. 



3 Again, note that M s (f) is denned for all s € Z>^, in particular for s > N, even if N is finite. 
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with H analytic and invertible near ttq, which implies that M s (£) is analytic near ttq if and only if di(7To) = 0. 
Finally, if N is finite, we have to make sure that M s (£) has no pole at a~ 1 iTN for s > N + 1. A necessary and 
sufficient condition for that is ^(cr ttiv) = 0. 

(c) Using the asymptotic conditions on the matrices m s (cr£) and m s +i(C), we obtain, as £ — > oo, 
M S (C) = m,K)2?(C)m^ 1 1 (C) 

I + O ^){ V v -^- k ){ XlC + o iC } \2( n + o(( n - 1 )) ( C o C° fe )( / + °@ 
Kl ' c" + o(c™~ 1 ). 



K 2 

Since M s (£) is entire by part (b), this completes the proof. □ 

3.3. To conclude this section, we remark that the method by which we have obtained the second equation 
of the Lax pair can be applied to derive a (second-order) difference equation for the orthogonal polynomials 



corresponding to the weight function ui. More precisely, one proves, by the same argument as Theorem 3.1, the 
following 



Proposition 3.3. Under the assumptions of Theorem 3. 1 (b ), the function M (£) = mx(o-£)D(£)m% (£) is entire 



The analogue of Theorem \3. J\ (c) also holds for the matrix M(Q. 

Now we have mx(o-() = M(£)m;a:(C)-D(C) -1 - Considering the (1, 1) and (2, 1) elements of both sides and using 
the explicit formula for mx (C) given in Theorem |2.4| , we obtain a system of equations of the form 

P fe «) = (M u (C)P fe (C)+cM 12 (C)P fc -i(C))di(C)- 1 , (3-5) 

C -P fe -i(0 = (M 2 \C)^(C)+cM 22 (C)P fe -i(C)H(Cr\ (3.6) 

where c = (Pfe-i, Pfc^i)" 1 and the M u (£) are the elements of the matrix M(Q. If we find c • Pfc_i(C) from the 
first equation and substitute it into the second one, we will get a relation between Pfc(C), Pfc(cC) and Pfc(cr 2 £). 
Even though the functions M lJ (£) may involve unknown parameters, one might hope that the equation we obtain 



in the end will only involve known parameters. We will work out an explicit example in subsection 11.2, for the 
weight function corresponding to the Charlier polynomials. We will see that the relation we get is exactly the 
standard difference equation for Charlier polynomials. 

4. Recurrence relation for Fredholm determinants 

4.1. We remain in the general setting of §H[ Thus we consider a discrete locally finite subset X = {it x } x=Q Q ^> 
where N £ Z>o U {oo}, and we are also given a strictly positive weight function lj : {x £ Z>o \x < N} — » C whose 
moments are finite. For all s £ Z> , s < N, we let 

3, = K}*=o, 2). = £\3.. 

Finally, we fix a natural number k < N and let 



m ^0- ^ m 21(£) m |?( C ) 

be the unique solution of the discrete Ricmann-Hilbcrt problem (X, w) with the asymptotics of diag(C fe , C~ fe ) a t 



infinity provided by Theorem 2.4 



Let a, {3 : {x £ Z>o|x < N} — > C be two functions such that a{x)(3(x) = ui(x) for all x. Consider the following 
kernel on X x X: 

K{K X ,-K y ) = \ 7T X 7Ty (4.1) 

a(x)/3(x)(<j>'(ir x )ip(ir x ) - i>' {-k x )4>{-k x )) , x = y, 
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where 0(C) = m£(C) = ffc(C) and V(C) = ro| x (C) = (Pk-i, Pk-i)^ 1 • flb-l(C)- Up to the factor of a(x)0(y), this 
is precisely the Christoffel-Darboux kernel for the family of orthogonal polynomials corresponding to the weight 
function u>. We will not need the specific form of the functions a and (3 in our computations. Note also that 
changing a and (3 while keeping their product fixed results in conjugation of the kernel K and hence has no effect 
on the Fredholm determinants studied below. 

4.2. For each s G Z>fc, s < N, we let K s be the restriction of K to 2) s x 2) s . We also denote by K and K 8 the 
operators on l 2 (X) and Z 2 (2) s ) defined by the kernels K and K s , respectively. The main goal of our paper is to 
derive a recurrence relation for the Fredholm determinants 

D s = det(l - K s ), s e Z> fc , s < JV. 



The resulting equation will be in terms of the elements of the matrices A s and M s (£) from the Lax pair (3.3) 



(3.4). Note that the Fredholm determinants are always well defined because if is a finite rank operator, as follows 



from the Christoffel-Darboux formula (see, e.g., JL8[): 



a(x)f3(y) ~ ^ {P m ,P m )u 



For a probabilistic interpretation of these determinants, see subsection 4.5 below 



Lemma 4.1. For each s G Z>/ c , s < N , the operator 1 — K s is invertible, and D s = det(l — K s ) ^ 0. 

Proof. The operator K is the orthogonal projection onto the subspace of l 2 (X) spanned by the polynomials 
Po: ■ ■ • 5 Pk-i- In particular, it has finite rank, its only eigenvalues are and 1, and the eigenvectors corresponding 
to the eigenvalue 1 are linear combinations of Pq, . . . ,-Pfe-i) i.e., the polynomials of degree < k — 1. It follows 
that the operator K s also has finite rank, hence 1 — K s is not invertible if and only if 1 is an eigenvalue of K s . 
Suppose that k < s < N and there exists / G l 2 (SQ s ) such that K s f = f. We extend / by zero to X, and we 
denote the extension by / G l 2 (X). Now Kf\^ = K s f = f = /| and =0. This implies that Kf\^ = 

because otherwise we would have ||if/|| w > \\f\\u, which is impossible because if is a projection operator. This 
shows that Kf — /, whence by the remarks above, / is a polynomial of degree < k — 1. But it vanishes on the set 
3s of cardinality > k, which implies that / is identically equal to zero. Hence, 1 — K s is an invertible operator, 
and det(l - K s ) ^0. □ 

4.3. For each s G Z>fc, s < TV, we have the unique solution m s (C) of the DRHP (3s,w|j ) having the same 
asymptotics at infinity as mx{Q- As in §[|, we assume that there exist entire functions di(Q, ^(C) such that 

di(7r )=0 and W ( X ~ ^ = n ■ for aUl < x < N. (4.2) 

uj{x) d 2 (Tr x ) 

The assumption that d2{o~~ 1 'K^) = if N is finite is not essential for us now since it was only used in §|3| to show 
that the function M S {Q is entire for s > N, whereas here we are only interested in the case k < s < N (cf. 



Remark |3.2|). Note that d 1 (n x )d 2 (Tr x ) ^ for all 1 < x < N. We let 

^(C)=(^ C) d2 ( C) ) and M s (0 = m s (o-OD(Om^ 1 ((). 



By part (b) of Theorem |3.l| , the function M S (C) is entire for all k < s < N. Using part (a) of this theorem, we 
obtain the following general form of the Lax pair for the solutions m s (C): 



»u,iO [l+-Ar)m.(Q, U ( ) . ,r (4.8) 



m.K)=M,(C)m a+ i(OI?(C)- 1 . 



(4.4) 
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In §[5], we will explain how one can obtain a recurrence relation for the matrix elements of M s (£) and A s . In the 
present section, we assume that the function M S (Q is known, and derive a recurrence relation for the Fredholm 
determinants using this function and the parameters p s , q s . 
For all s, k < s < N, we put 



m s (TT s 



M s (, s+1 ) = ( % % ) , M^ s+1 ) = ( % 



,22 



(4.5) 



where M' s (() = -^M s ((). Then we have the following recurrence relations for the matrix elements m* 1 and for 
the Fredholm determinants D s in terms of the parameters p s , q s , /if, v l s 3 . 

Theorem 4.2 (Recurrence relation for Fredholm determinants). Assume that p s ^ for all s > k, and hence 
also q s , r s 7^ for s > k. 

(a) For each s £ 7L>k, s < N, we have 

(b) For each s 6 Z>fc, s < N — 2, we have 

D s+2 D s+ i _ lu(s) ■ u s ■ (m] 1 ) 2 



11 



(4.6) 



D 



J s+2 



r] ■ di(x s+1 ) ■ d 2 (ir s+1 ) '' 



(4.7) 



wher 



Qs Qs 



for all s. 

Remark 4.3. We will show later (see Proposition [T| and Remark j6.3|) that under the assumptions of subsection 



1.2, p s is nonzero for all s > k. 



Proof, (a) First we take the residues of both sides of (4.2) at £ = ir s and use the jump condition on the LHS: 

lim m s+ i(C)w(7r s ) = Asm s (TTs)- 

Since the first column of the matrix w(ir s ) is zero, we deduce that the first column of the matrix ^4 s m s (7r s ) is 
zero, whence m 21 = — (p s /Qs)ml ■ Next we substitute £ = tt s +i in (4.4) and rewrite it as follows: 

m s+ i(7r s+ i) = Af7 1 (7r s+ i)m s (7r s )Z)(7r s+ i). 

Since detM s (() = detD(() by Lemma 2.1, the last equation can be written explicitly as 



,11 



,12 



■s+1 

91 22 

ro&i mf +1 



s+1 _ 



"Ms 1 f-l 1 



di(7r s+ i)d 2 (7r s+ i) 
Equating the (1,1) and (2,1) elements of both sides yields 



21 29 

mi mi z 



di(-K s+1 ) 
d 2 (n s +i) 



rn 



j / \— 1/ 22 11 12 21\ 1 / \-i / 22 , Ps 12 

V Qs 



,11 



and 



^ + "Ms 1 
Qs 



(b) For all s £ Z>fc, s < N, we put i? s = if s (l — ICs) 1 (note that this is well defined by Lemma LI), so that 
R s + l = (l -K a )~ l , and hence [R s + l)(s,s) = det(l - K s+1 )/det(l - K s ), i.e., 

Ds+i 



R s (s,s) 



Ds 



1. 



(4.8) 
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Now Theorem 2.3(h) in 0] gives (Situation 2.2 ibid, explains why this theorem is applicable here) 

R s (s, s) = g t (s)m~ 1 (ir s )m' s (ir s )f(s), 

where f(s) = (a(s),0)' and g(s) = (0, — /3(s))'. For the purpose of our calculation, we rewrite this formula as 
follows: 



R s (s, s) = -w(s) ■ ejm, 1 (7r s )m' s (7r s )ei, 
where &\ = (1, 0)' and e 2 = (0, 1)*. Similarly, we get 

R s +i(s + 1,8 + 1) = -w(s + 1) • 4™r+i(^*+0 m *+i( 7r *+i)ei- 



We substitute £ = tt s +i into (4.4) again: 

m a {ir s ) = M s {ir s+1 )m s+1 (ir s+1 )D(n s+1 )~ 1 , 



and we differentiate (4.4) at ( = ir s+ i to obtain 

d 



m' s (n s ) = t) 1 



(M^Om.+iCOPCC)- 1 ] 



(4.9) 
(4.10) 
(4.11) 

(4.12) 



If the derivative in (4.12) falls onto the third factor, the contribution of the corresponding term to the RHS of 
© is (by fill )) 

-rr^(s) • 4£>(7r s+ i)^| f= ^ +i OD(C)- 1 ) e x = const • (0, 1) ( * j)fjl=0, 

If the derivative in (4.12) falls onto the second factor, the contribution of the corresponding term to the RHS of 
© is, by ( PH ), Q and ( p0| ), 

-?7 _1 w(s) • e2L'(7r s+ i)m7_j 1 1 (7r s+ i)m / s+1 (7r s+ i)D(7r s+ i)" 1 ei 

-1 / v^CtTs+i) t -1 / \ / / \ 

= -V W ( S )t7 T ' e 2 m s+l( 7r s+l) 7n s+1 (7T s +i)ei 

= -w(s + 1) • 4™^(^+l) m s+l(>VK) e l 

= R s+1 (s + l,s + l). 

Therefore, if we subtract R s +i(s + 1, s + 1) from both sides of fl4.9| ), we get the following equation: 
R s+ l(s + 1, s + 1) - # s (s, s) = ■q~ 1 uj(s) ■ elm~ 1 (ir s )M' s (Tr s+1 )m s+1 ('ir s+1 )D(ir s+l y 1 ei. 
Since det m s = 1 for all s by Lemma |2.l| , the RHS of ( 4.13| ) equals 

mlXx mil, \ ( d x {-K s+x )- 1 



(4.13) 



fMs)'(o,i) 



,22 



,12 



"s+1 '"s+l 



d 2 (7T S +l)^ 



= ^-^(^di^+i)- 1 ^ 21 ™^^^! + ,fm>^ - ffmfm^ - v] 2 mfmf +1 ). 

Substituti ng th e formulas for mjj 1 , m ^-, , OTg^ derived in part (a) into the last expression and using (4.8), we 
hnd that (|4.13|) is equivalent to □ 



Corollary 4.4. FFif/i i/ie notation of Theorem {.i, we have 

D s+3 D s+2 \ _ (/if + (p s /q s ) ■ til 2 )' 



D s+2 D s+ i J r\ ■ di(7r s+2 ) • ^2(^+2) 
/or aZZ s 6 Z>fc, s < N — 3. 



• u s +i 



s+2 



D 



s+l 



D 



s+l 



(4.14) 



Proof. Immediate from (4/7) and (4.6). 



□ 
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4.4. It turns out that it is possible to extend the results of Theorem [4J2] to some cases where the orthogonality 
set X is discrete but not locally finite, i.e., has an accumulation point. The motivation for such an extension is 
the following. Even though the DRHP (X,w), as it is stated in §0, is not wen posed if X is not locally finite 



(for then we need to impose additional conditions near the accumulation point of the poles), the definition (4.1) 
of the kernel K still makes sense: the polynomials (f>(() = Pfe(C) an d VKC) = (Pk-i, -Pfc-i)^ 1 ' -Pfc-i(C) are well- 
defined. In fact, there exist several classical families of basic hypergeometric orthogonal polynomials for which 
the orthogonality set is discrete but not locally finite (see e.g. |l3|, Chapter 3). On the other hand, the solutions 
m s (C) of "restricted" DRHPs, and hence all the quantities derived from them, are also defined for a non-locally 
finite orthogonality set X, because their definitions involve only finite subsets of X. In particular, we can still 
consider the corresponding Lax pair as in subsection 4.3 and the scalar sequences p s , q s , m l s J 



(4.3) and (p 



Hi 3 , v l s 3 defined by 



It is therefore natural to ask whether the recurrence relation (4.7) remains valid in the case 



where X is not locally finite. We will see in subsection 4.6 that it does 



4.5. Let us recall the following probability-theoretic interpretation of the Fredholm determinants D s , see e.g. 
p9[ . In general, if X is a discrete, not necessarily locally finite subset of R of cardinality N + 1 (N 6 Z>o U {oo}), 
u> : X — > R is a strictly positive weight function whose moments are finite, {-Pn(C)}^=o ^ s the corresponding family 
of orthogonal polynomials and fc is a natural number, fc < N, we can consider a probability distribution *}3 on 
the set of all subsets of X of cardinality fc, defined by 



*P({ 



Xi, . 



I [ (Xi-Xjf- J[u)(Xi). 



(4.15) 



l<i<j<fc i=l 

Here Z is the unique constant for which the measure of the set of all subsets of X of cardinality fc is equal to 1 : 

k 

Z = E II (xi-x^-Hojixi), (4.16) 

{ai,...,x fc }C£l<i<.7<fc i-1 



Now if K denotes the kernel (4.1) defined in subsection [4.1^ then for any subset 2J C X, we have 



,Xk}) 



det(l = £ 

mm 

(the sum on the RHS is taken over all subsets {x\, . . . , Xk} C X of cardinality fc which are disjoint from 2)) 



(4.17) 



4.6. We now assume that X = {tt^I^q C R is discrete but not necessarily locally finite. We consider, as before, 
the subsets 2) s = {Tr x }%L s of X, and we are interested in the sequence {D s }<^ k defined by 



D s = det(l - K s ), K S = K\ 



2)»x2) s 



where K is the kernel ( pLip . The proo f of Lemma [4.1| is still valid and gives D s ^ for all s £ 1>>k- Thus we see 
from the discussion of subsection 4.4 that both sides of (4.7) are at least well defined in this situation. 



Proposition 4.5. All formulas of Theorem J^.i. and Corollary (.4 remain valid in the present situation 



Proof. It is o bvio us from the proof of Theorem 4.2 that (| 4.6h rem ains valid. Now let us fix s 6 Z>fc and prove 
that formula (4/7) also holds in the present situation (then ( |4. 14|) follows automatically). Let L £ Z, L > s + 2. 
We write X^ = {tt x }£ =0 C X. Since is finite, all of the discussion of subsections [I. is valid for X (L ^ in 
place of X. Since L > s + 2, it is clear that replacing X by X^ (and keeping the same weight function) has no 
effect on m s (C), p s , q s , m l s J ', jU* J ', v\? ' . But of course, the quantities D s do change. Let denote the Fredholm 

determinants defined as in subsection 4.2 for X^ in place of X. Then Theorem 4.2(b) gives 



D 



(L) 
s+2 



D 



(L) 
s+1 



u)(s) ■ u s ■ (ml 1 ) 2 

77 • dl(7T s+ i) • d 2 (7Ts+l)' 



D 



(L) 
s+1 



D 
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It remains to observe that D { s+1 /Ds ' — D s+ i/D s for all s. Indeed, the discussion of subsection L5 gives the 
formula 



0<ii <i2 < - <ifc<s — 1 



(4.18) 



where is given by ( 4.15| ). If ^( L ) denotes the probability distribution on the set of all subsets of X^ of 
cardinality k defined similarly to ^p, then we also have 



0<h<i 2 <— <i fe <s-l 



(4.19) 



We note that the summations in ( l.lq ) and (4.19) are over the same index set, and the only difference between 
the definitions of and <p( L - ) is in the normalization constant Z. Of course, when we take the ratios -DS/ 
and D s+ i/D s , the normalization constants cancel each other, completing the proof. □ 

5. Compatibility conditions for Lax pairs 
5.1. In this section we study the compatibility conditions for the Lax pairs of the form considered in §|| §[|: 



l-s+l 



(C) = U + 



A, 



c 



,(0. 



(5.1) 



m s «) = M a (0ro a+ i(C).D(Cr 1 . (5-2) 

The general notation and conventions are those of §[$| §[|. As in the second part of we do not assume that the 
orthogonality set X is locally finite: we have already remarked in subsection 4.4 that all of our arguments related 
to Lax pairs only involve finite subsets of X. 

Lemma 5.1. Fix s £ Z> , s < N, let A 6 Mat(2,C), and define m(£) = (J + (C - n s )- 1 A)m s ((). Then 
m(C) = ms+l(C) an d on ly if the matrix A satisfies the following two conditions 



A ■ m s (-K s )w(iT s ) = 0, 
m s (7r s )w(7r s ) + A ■ m' s (ir s )w(Tr s ) = A ■ m s (n s ). 

In particular, for a fixed s, there is a unique matrix A satisfying (pTjj and ft5^2), namely, A = A s 



(5.3) 
(5.4) 



Proof. By the uniqueness of m s +i(C), we only have to verify that m(£) satisfies the same residue conditions as 
m s +i(C) if and only if (5.3) and ( |5.4| ) hold (note that the asymptotics of m{Q and m s +i(C) as C —* oo are clearly 
the same). Now ifO<x<s — 1, then since (/ + (C — 7r s ) _1 A) is analytic near ir x , it is clear that the residue 
condition at ir x for m s (Q implies one for m(Q. Thus, we only need to consider the residue condition at the pole 
C = 7r s . Since m s (Q is analytic near ir s , we have Res^ =Tr ,w(C) = A ■ m s (jr s ). On the other hand, the limit 

lim m(()w(ir s ) = lim (/ + (C - ■n s )~ 1 A)m s (Qw('n s ) 

exists if and only if ( |5.3j ) holds. Moreover, if ( |5~3| ) holds, this limit equals m s (w s )w(n s ) + A ■ m' s (n s )w(n s ), so 
m(£) satisfies the required residue condition at £ = w s if and only if (5.4) holds. □ 



Theorem 5.2 (Compatibility conditions for Lax pairs). Fix s G 7L>k, s < N — 1. 
(a) We have 



A s 



M S (C)=M S+1 (C) / + 



.4 



s+1 



(5.5) 
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(b) Conversely, assume that M : C 
and 



Mat{2, C) is an analytic function, A € Mat(2, C) is a nilpotent matrix, 



M a (C)=M(C) 

11). ThenM(Q 



A 



(5.6) 



A 



s+l- 



C - 7Ts+l . 

M s+ i(C) and A : 

D- 

Remark 5.3. This theorem provides a recipe for computing A s+1 and M 8 +i(£) if we know A s and Jlf,((). Indeed, 
one simply needs to find the unique solution of the compatibility condition which satisfies A 2 S+1 = 0. 



where M a (£) and A s are defined by (5.1 ), 
Equation (|5.5|) is the compatibility condition for the Lax pair (5.1), (E 



Proof, (a) If we replace £ by <t£ in (5.1) and then substitute (5.2) into the result, we obtain 

A. 



m s+1 (a() = [I - 
On the other hand, if we substitute (fO) into ( 



m s+ i(aC)=M s+1 (C) 1 + 



M s (0m s+ i(QD(0- x . 

1, we get 



and then replace s by s 
A s+ i 



s+l 



(5.7) 



(5.8) 



Comparing (|5.7| ) and (|5 
(b) From (O) and (O), 



) yields fl5_5|). 
we find that 



M(C) / + 



.4 



C 



= M S+1 (C) / 



,4. 



s+l 



Using Lemma 2.1 and ( |5.2| ) with s replaced by s - 
invertible near ir s+ i. Now since A 2 = 0, we can rewrite (5.9) as 



C - 7T S +1 

we see that detM s+ i(£) = 



A 



s+l 



A 



C - 



(5.9) 

detD(C), and hence M s+ i(() is 

(5.10) 



The RHS is analytic for £ 7^ 7r s +i, and the LHS is analytic near 7r«+i. Thus, both sides of ( 5.10| ) are entire 
functions. But the RHS tends to I as £ — > 00, so by Liouville's theorem, both sides are equal to / for all (. This 
proves that M(Q = M s +i(C) and A = A s+V □ 

5.2. Recall that the formulas of Theorem [4.2j have been derived under the assumption that the parameter p s 
does not vanish for all s > k. Let us now establish the no n- va nishing of p s for the weight functions lo and 
orthogonality sets X satisfying the assumptions of subsection L2. We need the following well-known fact. For 
the reader's convenience, we also provide a proof. 

Lemma 5.4 (Zeroes of discrete orthogonal polynomials). Let 3 C C be a finite subset, and assume that 3 is 
contained in a closed interval [a,b] C R C C such that a is the minimal element of 3 and b is the maximal 
element of 3- Let ui : 3 —> K oe a strictly positive weight function. Then there exists a unique family {P n (()}n=o 
of orthogonal polynomials corresponding to to, where L = card(3) — 1. The coefficients of each P n (C) are real. 
Moreover, for any < n < L, all zeroes of P n (Q) are real and are contained in the open interval (a, 6). 

i<d 



Proof. Since u) is strictly positive, the restriction of the corresponding inner product (•, -) u to the space 
of real polynomials of degree at most d is nondegenerate for each < d < L. Thus, there exists a unique family 
of real orthogonal polynomials {P n (C)}n=a corresponding to ui, and {P n ,Pn)u 7^ for < n < L. Now fix 
1 < n < L, and assume that P n (() has fewer than n zeroes in the open interval (a, b). Let z\,...,z m be the 
zeroes of P n (Q in (a, b), listed with their multiplicities, and let Q(Q = (C — z i) 1 • • (C — z m)- Then, since P n (() 
and Q(C) are real polynomials, we have either P n (C)Q(Q — f° r ah C G [ a -, b], or P n (QQ{Q < f° r all C G [ a i ki- 
ln addition, since the degree of P n {C) is l ess than the cardinality of 3, there exists z G 3 such that P n (z)Q(z) ^ 0. 
This implies that P n (C) an d Q(C) are not orthogonal with respect to to. Since the degree of Q(C) is less than that 
of P n {(), we have a contradiction with the definition of orthogonal polynomials. □ 
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Now we can prove 



Proposition 5.5. 

7T > 7Ti > 7T 2 > • • 



With the notation and conventions of 
. Then p s =/= for all s > k, s < N . 



and 



assume that either ttq < tv\ < tt2 ■ 



or 



Proof. Fix s > k, s < N. By Lemma 5.4, there exists a family {P n {C)}n=o °f polynomials orthogonal on 
{ttq, . . . , 7r s _i} with the weight function given by the restriction of u to {0, 1, . . . , s — 1}. Moreover, as 7r s lies 
outside the interval between ttq and tt s ~i, we have P n (jT s ) ^ for all < n < s — 1. Now since k < s — 1, we 
know from Theorem 2.4 that the first column of the matrix 77i s (£) has the form (Pfc(£), c • Pfc_i(C))\ where c is 



a nonzero constant. On the other hand, by Lemma 5.1, we have 



m s (Tr s )w(ir s ) = A s ■ (m s (7r s ) - m' s (7r s )u!(7r s )). 



(5.11) 



If p s = 0, then because the matrix A s is nilpotcnt by Theorem 3.1(a), we have either q s = or r s — 0, i.e., either 
the first or the second row of A s is zero. By ( 5.1l[ ), this implies that either the (1, 2) or the (2, 2) element of the 
matrix m s (Tr s )w(Tr s ) is zero. This contradicts Pfe(7r s ), p c _i(7r s ) =/= 0. □ 



6. Initial conditions for recurrence relations 



6.1. In this section we derive the initial conditions for the recurrence relations (4/7) and ( |5.5| ). We keep the 
general notation and conventions of §|| and §0. Recall in particular that k is a natural number that controls the 
asymptotics at infinity of the solutions of all DRHPs that we consider. As in we do not assume that the 
orthogonality set X is locally finite. 



Proposition 6.1. The solution mfe(C) of the DRHP ({no, 
m k{C) ~ ( C c °k ) as C -> oo is ^wen &?/ 

(C _ 7ro ) (C _ 7ri )...( C _ 7rfe _ 1 ) 



,7T/c-l}, W 



mki0 ~ [ (C - TToXC " TTl) • • • (C " 7T fe _!) E™=0 ^ 



{7ro,...,7r fc _i}' 







wii/i ^/ie asymptotics 



(C-TrorHC-Tri)- 1 -"^-^-!)- 1 



where 



for all < m < fc — 1 . 



: w ( m ) 1 ■ n ^ 



(6.1) 



(6.2) 



0<j<k-l 



Proof. Let m(C) be the matrix defined by the RHS of (|Q|). It is clear that m(£) has the required asymptotics at 
infinity. Hence we only have to show that (6.2) is the (unique) choice of constants p m which makes m(Q satisfy 
the required residue conditions. Now if < x < k — 1, then the (2, 2) element of the matrix Res^^m^) equals 



n (' 



0<l<k-l 



and the (2,2) element of the matrix lim^^ Trx m(£)w(ir x ) equals 



0<l<k-l 



The other elements of both matrices are zero. Equating the last two expressions yields (6.2) 



□ 
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6.2. Now we use Lemma 5.1 to find the matrix A k . 



Proposition 6.2. The elements of the matrix 



Pk qk 
Tk ~Pk 



are given by the following formulas: 



( k-l 
qk = \Pk+Yl 



pm 



m=0 



where 



E 



fe-i 



0<j<k 



(6.3) 



p k = uj(k) 1 • Y[{TT k - 7Tj ) 2 ; 

3=0 



Pk = -qk ■ 



fe-i 

E 

rn— 



7Tfc — 7T„ 



and 



fe-i 



-«* • i E 



v m=0 



7I"fe - TTr, 



(6.4) 



(6.5) 



Remark 6.3. It follows from ( |6.2j ), (3.3) and ( |6.4| ) that if the orthogonality set X is contained in R and either 
7To > 7Ti > 7T2 > • ■ ■ or tto < 7Ti < 7T2 < ■ • ■ (and the weight function u> is strictly positive), then p m > for 
< m < fc — 1, qk > 0, and hence p k ^ 0. 



Proof. It follows from Lemma 5.1 that A k is the unique matrix satisfying the following system of equations: 

A k ■ m k (-K k )w{-K k ) = 0, (6.6) 
m k (-K k )w{-K k ) + A k ■ m' k (-K k )w{iT k ) = A k ■ m k (ir k ). (6.7) 



Substituting (3.1) into ( p.6[ ) yields (|6.4| ) . It remains to prove (6.3), for (6.5) then follows from p\ = —q k r k . 
To this end, we consider the (1,2) elements of both sides of (6.7). The first summand on the LHS contributes 



uj{k) ■ P3^ = q (7Tfc — nj) to the (1, 2) element. Now we consider the second summand. We can rewrite it as 

d 



d_ 

<K 



(A k • m fe (C)w(7r fc )) = — < (C - 7r ) • • • (C - TTfc-i) • w(fc) • A k 



1 



u Z-im=0 C-ir m / J 

If the derivative falls onto the factor (£ — ttq) ■■ ■ (£ — 7Tfc_i), the corresponding term is zero because of ( |6.6[ ). 
Hence the whole expression equals 

« — i n Y^k— i Pm 

U —qk ■ Z^ m =0 (7r fc -Tr m ) 2 



(7T fe - 7T ) ■ • ■ (7T fe - 7Tfc_l) ■ A k ■ Uj(k) 











U Z-tm = Q (TT k -TT 







Finally, the (1,2) element of A k ■ m k (-K k ) equals q k ■ 11^=0 (^fc — n j) 1 - Thus, comparing the (1,2) elements of 
both sides of (6.7) yields 



fe-i 

Yl^k-TTj) 

J=0 



k-l 



;( fc ) • E 



Pn 



k-l 



k-l 
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which gives (3.3). 



□ 



Remark 6.4. The two propositions we have just proved give explicit formulas for the matrices rrik(C) an d A k . 
Usin g these formulas, we can also find the functions mk+i(Q and M/ S (C)- Indeed, from the Lax pair (3.3) and 
(3.4), we have 

A k 



I 



C — TTfe 



wfc(C), 



M fe (C) = mkiaODiQm^it) = m k (a()D(t) m - 1 (0 / 



.4, 



(6.8) 



(6.9) 



Even though this gives an explicit formula for M k (C), it is cumbersome to use it in practice. In the case where 
the matrix Z?(£) is linear in £, a more explicit version is available: 

Proposition 6.5. Assume that rfi(0 = Ai£ + /xi, (^(C) = A2C + M2J where Ai, A2, /ii, /i2 G C are constants (some 
of which could be zero). Then 

?7 fe (AiC + Mi) +?7 fe Ai(7r -7T fc -jo fe ) -^AiOfe 
-r/~ fe A 2 r fc + fa*" 1 *! - rr fe A 2 ) £^ = X Pm T fe (A 2 C + M2) + rT fe A 2 (;Pfc + ^ fe - 7r ) 



M fc (C) 



(6.10) 



where p k ,q kl r kl p m are given by (6.4), (6.8), ^6.3j) and (6.i.) 



Proof. Since Mfc(0 is an entire function, it suffices by Liouville's theorem to show that (6.10) holds up to terms 
of order C^ 1 as £ — * 00. To that end, note that by (3.1) and (3.2), we have 

»? fc (C-Ti)--.(C-7r fc ) 



m k (a() = 



v ^-i( C _ Wl) ...( C _^)^ 0? _ +j 
But it follows from (p^9) that 



(6.11) 



M fe (0 = ^m fe «) 



C" fc 
c fc 



£>(0" WO 



C" fc 
c fc 



I- 



Ai 



Substituting (3.1) and ( 6.11 ) into the last formula, we arrive at ( 6.10| ). 



□ 



6.3. The final result of this section is the computation of the Fredholm determinants D k and D k +i- We use 
the probability-theoretic interpretation of the Fredholm determinants D s given in subsection 4.5. One can show 
(this is a standard random matrix theory argument, see e.g. [[l4|) that the constant Z given by (4.16) is equal to 
the product of the norms squared of the first k monic orthogonal polynomials: 

fe-i 

Z=Y[(Pi,Pi)u. (6.12) 

i=0 

Now we prove 

Proposition 6.6. With the notation of §0, let X C M be a discrete set, let {P n {()} be the family of orthogonal 
polynomials corresponding to a strictly positive weight function lo : X — > R, and let Z be given by (6.1i ). Then 

fe-i 



D k = 



\- n fa-^) 2 -ri<"(o. 



0<i<j<fe-l 

D k+1 = uu{k) ■ ql 1 ■ D k ■ J](7T fe - tt ; ) 2 , 



fc-1 



(6.13) 
(6.14) 



where q k is given by ( 6.S ). 
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Proof. Recall that for all s e Z>fc, s < N, we have defined 3s = {^I^O' — 3E\3«i an d D s — det(l— K\ x ^ ). 
Hence a subset of X is disjoin from 2) s if and only if it is co ntain ed in 3s - Th ere exists only one subset of 3fe of 
cardinality k, namely, 3fc = {tto, ■ ■ ■ , ^fc-i} itself. Applying ( [4.17 ) yields ( |6.13| ). Next, there are k + 1 subsets of 



3a-+i of cardinality k, namely, those of the form 3fc+i \ {"m} for < m < fc. Applying ( 4.17 ) gives 



2=0 



1 



u;(m) 



0<i<j<k 



Using @, (fj) and we see that the last equation is equivalent to ( |6.14| ). 



□ 



7. Lax pairs for discrete orthogonal polynomials of the Askey scheme 

7.1. In this section we specialize to weight functions appearing in the orthogonality relations for the hyperge- 
ometric orthogonal polynomials and the basic hypergeometric orthogonal polynomials of the Askey scheme. We 
use |D| as our main reference for the orthogonal polynomials of the Askey scheme. We are only interested in 
those families for which the orthogonality set is discrete. Since our ultimate goal is to derive a recurrence relation 
for the associated Fredholm determinants, we do not impose the local fmitcness condition on X. However, the 
basic assumptions of subsection 1.2 have to be satisfied in order to use our approach (in its present form). These 
assumptions are not satisfied for the following families of discrete orthogonal polynomials: the Racah polynomials 
(p3[. §1.2), the dual Hahn polynomials (p3y, §1-6), the q-Racah polynomials (|L3|, §3.2), the big g-Jacobi poly- 
nomials (pq|, §3.5), the big g-Legendre polynomials ([^3|, §3.5.1), the dual g-Hahn polynomials (|R|, §3.7), the 
big g-Laguerre polynomials ([[l3], §3.11), the dual g-Krawtchouk polynomials (|l3), §3.17), the Al-Salam-Carlitz 
I polynomials (p3[, §3.24), and the discrete g-Hermite I polynomials ((ll|, §3.28). 

7.2. Now we list the families of hypergeometric and basic hypergeometric orthogonal polynomials for which our 
results do apply. Instead of writing out the whole Lax pair in each case, we only give the orthogonality set X, the 
weight function u>(x), the affine transformation a : C — > C, and the corresponding entire functions <ii(C), £fe(C) 
which satisfy the assumption of Theorem |3.l| (b). For the basic hypergeometric polynomials, we assume from now 
on that < q < 1. This restriction ensures that X C K and the weight function is strictly positive and has finite 
moments. 

• Hahn polynomials (||, §1.5): X = {0, . . . , N}, N e Z> ; 

, , fa + x\ f/3 + N - x\ 
uj{x) =1 I I I, where a, p > — 1 or a, p < —A; 



trC = C-l, d 1 (0=at-P-N-l), d 2 (t) = (t-N-l)({ + a). 
• Meixner polynomials (M], §1-9): X = Z>o; 



c x , where (3 > and < c < 1; 



ctC = C-1, di(0 = C, d 2 (0 = c( + c(/3 - 1) 
Krawtchouk polynomials (pf, §1.10): X = {0, . . . , N}, N e Z> ; 



N 



uj{ x ) = ( -' )p x (1-pY- x , where < p < 1; 



<rC = C-l, d 1 (Q = C, da(C) 



p-1 



(C-AT-1). 
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Charlier polynomials (M), §1-12): X = Z> ; 



a 

Lu(xj — — , where a > 0; 
or! 



aC = C-l, di(C) = C, d2(0 = o. 

g-Hahn polynomials (pi, §3.6): X = {g-^la; = 0, . . . , N}, N G Z> ; 



(ag;g) x (g N ;q) x , x _ x _jy 

w(x) = -— -— j^—iafiq) x , where < a,/3 < g or a,/3 > g JV ; 

{q;q)x[p q ,q)x 



aC = q(, d 1 (0 = af3((-l)(C-[3- 1 q- N - 1 ), d 2 (C) = (C - «)(C - 
Little g-Jacobi polynomials (JT^j, §3.12): X — {q x \x G Z> }; 

u(x) = <Khq]q } x (aq) x , where < a < g" 1 and b < g" 1 ; 

(<7;g)x 

<rC = « -1 C, di(C)=C-l, (fa(C)=a(6C-l). 
g-Meixner polynomials (jl3), §3.13): X = {g _2: |x G Z> }; 

w (a;) = -^-^ - c a: g(^, where < b < g" 1 and c> 0; 

[q;q)x{-ocq;q) x 

<j( = q(, di(C) = (C-l)(C + te). da(C)=c(C-6). 
Quantum g-Krawtchouk polynomials §3.14): X = {g _x |x = 0, . . . , AT}, TV G Z>o; 

/ \ ipq\q)N-x i txjv-s (Y) - — jv 

uux) = 7 r— 7 r (— 1) q K2 > , where p > q ; 

(q;q)x(q;q)N-x 

*C = ?C, rfi(C) = (C-i)(re A,+1 C-i), d 2 (C) = i-g Ar+1 C 

g-Krawtchouk polynomials §3.15): X = {g _a |x = 0, . . . , N}, N G Z> ; 

u( x ) = W , '? } (-p)- x , where p > 0; 

(<?;<?)* 

aC = ?C, di(C)=p(C-i), <fc(0=? _JV -«C. 

Affine g-Krawtchouk polynomials (|0|, §3.16): X = {g _x |x = 0, . . . , A^}, N G Z>o; 

/ \ (pq;q)x(q;q)N , , „ ^ ^ _ x 

= t — w — ^ — (w : where < p < q ; 
{q;q)x{q;q)N-x 

^C = <zC, di(C)=P(C-i), d2(C) = (C-p)(9 w+1 C-i)- 

Little g-Laguerre /Wall polynomials (g|, §3.20): X = {g*|x G Z> }; 

(ag) 2 

u(x) = t r— , where < a < g ; 

a( = q- 1 (, rfi(C) = C - 1, rf 2 (C) = -a- 
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Alternative g-Charlier polynomials (jL3|, §3.22): X — {q x \x G Z>q}; 

u(x) — 7 \ — Q 2 > i where a > 0; 

aC = g- x C, di(C) = C-l, rf 2 (C) = --C 

g-Charlier polynomials ((l3), §3.23): X = {g _x |x G Z> }; 

= 7 r-?^ 2 -', where a > 0; 

^C = <zC, di(0 = C-l, d2(0 = a- 
Al-Salam-Carlitz II polynomials ( |13]| , §3.25): X = ^ | a? G Z>o}; 

= — : — , where a > 0; 

(q;q) x {aq; q) x 

a( = q(, d x (0 = (C-l)(C-a), d 2 (C) = o. 
7.3. The next three sections (§§|b| fiol) deal with the various possible ways of "solving" the compatibility con- 



ditions for the Lax pairs listed above. By a "solution" of a compatibility condition of the form (^j) we mean 



a collection of formulas which allow us to express the entries of the matrices A s+ i as rational functions 

of the entries of the matrices A#, A s , where M S (C) = AfjV + • • • + Afj 0) , A# G Mat(2,C), for all s. The 
most general case where we have been able to solve the compatibility condition explicitly is the one where the 
functions di(() and ^(C) are either linear or constant; this is described in §|[ The resulting formulas can be used 
for practical computations, but they do not appear to be related to any known systems of difference equations. 
In certain more specialized cases we have been able to reduce the compatibility condition to one of the equations 
of H. Sakai's hierarchy in [ p7| . 



7.4. In §|9| we solve the compatibility condition (^g) in the case where the orthogonality set has the form 
X = {x G Z> |a; < N} (N G Z> U {oo}) and the functions rfi(C) and d 2 (C) are linear, by a method different 
from the one used in §|^ . We show that if both d\ and e?2 are nonconstant, then the compatibility condition 
is (generically) equivalent to the d — Py equation of H. Sakai jlffl, and if di is constant, then the compatibility 



condition is (generically) equivalent to the d—Piv equation ibid, (see Theorem ^J). This result allows us to write 
down explicit solutions for the recurrence relations corresponding to the Meixner polynomials, the Krawtchouk 
polynomials, and the Charlier polynomials — see § fll| . 



7.5. As for the basic hypergeometric orthogonal polynomials, we show in § |10| (see Theorem 10.1 (b)) that in 
the case where the orthogonality set has the form X = {(? _s }fLo (N €E Z>o U {oo}) and the functions di(C) 
and <i2 (C) are linear and nonconstant, the compatibility condition for the corresponding Lax pair is equivalent 
to the q — Pyi system of M. Jimbo and H. Sakai, for a certain choice of the parameters. Since the case where 
X = {q s \s G Z> } is reduced to the former one after replacing q by this situation suits the following 
families of basic hypergeometric orthogonal polynomials: the little q-Jacobi polynomials and the q-Krawtchouk 
polynomials (it does not suit the alternative g-Charlier polynomials because the constant terms of the functions 



and d2(C) must be nonzero, cf. Theorem |l0.3| ). If one of the functions di(C) and <fe(C) is linear and the 



other one is constant, it is possible to reduce the corresponding compatibility condition to a degeneration of the 



q — Pyi system. This process is described in subsection 10.3. It allows us to solve the compatibility conditions 
for the little g-Laguerre/Wall polynomials and the g-Charlier polynomials. 

It turns out, however, that the method of solving the compatibility condition by reducing it to the q — Pyi 
system (or its degeneration) is rather difficult to carry out in practice. So to find a recurrence relation for 
the Fredholm determinants associated to classical families of basic hypergeometric orthogonal polynomials, we 
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prefer to use the more general formulas of §|| (see §[Tl|) . The disadvantage of the formulas of §|[ as compared to 
q — Pyi, is the fact that the recurrence step substantially involves more than two sequences, while for q — Pyi 
two sequences suffice, cf. Theorems 8.2 and 10. 3| below. 



8. Solution of the compatibility condition: the general case 



8.1. In this section we "solve" (in the sense of subsection 7.3) the compatibility condition 



s+1 



M s (C)=M s+1 (0 (1 + 



-4, 



c - k s 



(8.1) 



derived in §|] (cf. equation ( |5.5| )), in the case where the matrix D(Q defined in Theorem |3.l| (b) depends linearly 
on £. Our method is based on the following simple observation. 

Proposition 8.1. If M s (£) = A • £ + C s for all s, where C s does not depend on Q and A is a fixed matrix 
independent both of £ and of s, then under the assumption that A 2 S+1 = 0, the compatibility condition (8.1) is 
equivalent to the following system of linear equations: 



(Tr s+1 A + C s +n [ A S A) • A s+1 = r, 1 A s ■ (tt s+1 A + C s ), 
C s+ i =C S + v^AsA - AA s+ i. 



(8.2) 
(8.3) 



Proof. C omp aring the asymptotics of both sides of (ST) as £ — » oo yields fl8.3| ). If we take the residues of both 
sides of ( |S.l| ) at £ = tt s +1) we obtain 



T) 1 A S ■ (tt s+ iA + C a ) = (tts+iA + C a+ i) ■ A s+1 . 
Substituting (^3) into the last equation and using the assumption that A 2 S+1 = gives (8.2). 



□ 



8.2. We now note that if the matrix 7r s +iA + C s + r\ A S A is invertibl e, th en the system fl8.2|) , (8.3) already 



has a unique solution (A s+ i,C s +i). Since the compatibility condition (8.1) always has a unique solution by 
Theorem 5.2(b), it follows from Proposition 3.1 that in this case, the solution of (3.2), (8.3) is also the solution 
of (3.1). Even though we cannot prove that the matrix 7r s+ iA + C S + rj~ 1 A s A is invertible in general, the explicit 
computations we have carried out for five families of basic hypergeometric orthogonal polynomials (see j |Tl] ) show 
that the result below (Theorem pT2) has practical significance. 



8.3. We introduce the following notation. Suppos e th at di (Q = AiC+Mij ^(C) = ^2C+/ i 2, w here Ai , fit , A2 , jig G 
C are constants. Then it follows from Theorem 3.1(c) that the assumption of Proposition 8.1 is satisfied for 
A = diag(«i, K2), where k± — rj k \i and K2 = r/ _fe A2. Let us write 



A s 



r s -p s 



and C s 



ot s f3 s 



Finally, define e s = det(7r s+ iA + C s + rj 1 A S A). Now we can state 
Theorem 8.2. We have 

e s = d 1 (TT s+1 )d 2 (TT s+1 ) + 77 _1 Ki(p s (5 s - r s [3 s ) - n^ 1 K 2 (Psa s + q s J s )- 



(8.4) 
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If e s 7^ 0, then the following formulas hold 
Ps+i = -ry"V7 le ; 

Qs+i 
r s +i 



{ r> 
■ (PsPs 


+ q s S s + K 2 ir s+1 q s ) ■ (r s a s - p s j s + Ki7T s+1 r s ), 


(8.5) 


= fl^ls 




■ (p s f3 s + q s S s + K 2 7T s+ iq s ) 2 , 


(8.6) 






(r s a s -p s "( s + KiTr s+ ir s ) 2 , 


(8.7) 




= Us - 


f rj^KiPs - Kip s+ i, 


(8.8) 


Ps+1 


= Ps- 




(8.9) 


Js+1 : 


= Is - 


f Tj^K\T s - K 2 r s+ x, 


(8.10) 


S a +l 


= 6 S - 


- if 1 n 2 p s + Kzps+i, 


(8.11) 






(kiS s - K 2 a s ) + -| • KiPs, 


(8.12) 




<ls 


Ts 



where u s is defined in Theorem \4.q (b). 

We omit the proof, as it consists entirely of straightforward computations. We first derive ( |S.4| ) using the 
identity detM s (£) = det-D(C) which follows from Lemma |2.l| . If e s ^ 0, we rewrite ( |3~2] ) as 

A s+1 = rf 1 • (tt s+ iA + C s + r^AsA)- 1 ■ A s ■ (ir s+1 A + C s ). 

Writi ng out the matrix product on the RHS explicitly yields (^-(jD!) ■ Then (|]| gives (^)-{ ^.ll| ). Finally, 
(3.12) follows immediately from the definition of u s . 

9. The fifth and the fourth discrete Painleve equations 

9.1. In this section we assume that the orthogonality set is of the form I={i6 Z>o|x < N} (N s Z>oU{oo}), 
so that, with the notation of §[| cr£ = C — 1 and rj = 1. Our goal is to prove that if the functions and cfe(C) 

are linear, then the compatibility condition (5.5) is equivalent to either the fifth or the fourth discrete Painleve 
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Recall from Theorem 3.1 that since ttq = 0, we must have d\(0) = 0, and because only the ratio di(C)/c?2(C) 
matters, we may assume without loss of generality that rfi(C) = C> and write cfe(C) — £C + T f° r some £, r e C 
(unless otherwise explicitly stated, we do not exclude the possibility £ = 0). By Theorem 3.1(c), we can write 

a 11 



M S (C)=AC + C S , where A 



1 

£ 



and C s 



C 



21 



Cf 

cf- 



Then the compatibility condition (5.5) takes the form 

A 



As 



C-(s + 



C-(s + l) 



(9.1) 



Yy) '(AC + C S ) = (AC + C S+1 ) • (i- 

where A s — f _ 9 p s J . The following result allows us to find a convenient reparameterization of the matrices A s 
and C s which leads to an explicit solution of (9.1). 

Lemma 9.1. We have 

C 1 s 1 +ps = -k, Cf-£ Ps =£fc + r, CfCf = Cl l Cf (9.2) 

for all s £ 1>k, s < N. 



Proof. Taking the asymptotics of both sides of ( ) as £ — > oo gives 

C s + A s A = Cs+i + AA s+ i, 

which implies that 

C] 1 + Ps = ClX, + p s+1 and Cf - &s = Cf +1 - £ Ps+1 



(9.3) 
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for all s £ Z>fc, s < N. This means that the expressions + p s and C s — £p s do not depend on s. But we 
know from Proposition 6.5 that 

C k = 



-k 



Pk 



-'Ik 



-£r k + (1 - E^Jo Pm r + £( Pfc + k) 



(9.4) 



whence C^ 1 



Pk 



-k and Cjj: — S,p k = £fc + r, proving the first two equalities in (|9.2| ). To prove the third 

□ 



equality, note that det M s (() — det -D(C) for all £, by ( J3.4D and Lemma |2.l| , and take £ — 

9.2. It follows from Lemma |9.l| that for all s £ Z>fc, s < N, the matrices A s and C s can be naturally parame- 
terized as follows: 



A, 



(k 



-1 -a s (3 s 
l/{a s p s ) 1 



C„ = 



(r - £& s )/& 



(9.5) 



Remark 9.2. This parameterization is taken from Q (see equation (6.8) ibid.). The proof of the theorem below 
is based on the same idea as the proof of Proposition 6.3 in Q|. In fact, the situation considered in §6 of [Q 
corresponds, with minor modifications, to the weight function for Meixner polynomials (see §0 above). The only 
essential difference with the present paper is that we consider a slightly more general situation by letting c?2(C) 
be an arbitrary linear function, which allows us to treat the cases of the Krawtchouk and Charlier polynomials, 
as well as the case of the Meixner polyn omi als. 

Note also that the parameterization (9.5) is only valid if the matrices A s and C s are su fficie ntly generic. If, 
for instance, C] 2 = 0, but C^ 1 / 0, then (3.5) does not make sense. When we try to solve (9.1) in terms of the 
parameterization (9.5), we will encounter a similar difficulty: the formulas will involve rational functions of the 
parameters, and it is not clear a priori that the denominators of all fractions do not vanish. We refer the reader 
to pj, §6, where this problem is discussed in detail. The argument of M can be easily adapted to our situation. 



Theorem 9.3. (a) Assume that £ ^ 0. Introduce new variables f s , g s by 

fs = -k-b s + - , g s 

I— a. 



Then with the parameterization (9.5), the recurrence relation (9.1) has the following solution: 

s t/£ + s + 1 



ft 



+i 



u 



s+l 



9s+i9s 



-(fe+r/o + T— + / : ' 7 

1 + g s 1 + t,g s 
(f s+1 -l-s)(f s+1 -l-s + k) 

Zf s+ l(f s+ l+k+T/0 

(1 + g s+1 )f s+1 + (k + t/Ci9s+\ - s - 1 



Ps g s +i (l + £# s )/ s+ i 

(b) Now let £ = 0, and introduce new variables f s , g s by 



fs 



ra s 



1. 



Then with the parameterization (9.5), the recurrence relation (9.1) has the following solution: 

rg s 



fsfs+l — 



{9s 



9 s+i = 



0s 



s+l 

fs + 1 
T 



1 



1) 



k + 2s + 3. 



fs{g s 



1)' 



(9.6) 

(9.7) 
(9.8) 
(9.9) 

(9.10) 

(9.11) 
(9.12) 
(9.13) 



:2l i 
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Remark 9.4. (a) If we set / = / s , / = f s +i, 9 = g s , 9 = 9s+i, then the relations (9.7), ( |9.8| ) form a special case 
of the difference Painleve V equation (d — Py) of JT7j ], §7. The parameters A, ao, a±, 02, &3, &4 in our case are 
as follows: 



a o = r /£ + s + 1, ai 



a 2 



s, 



a 3 = — (fe + t/£), a 4 = k, A = ai + 2a 2 + a 3 + a 4 + a = 1. 

(b) If we set / = / s +i, f — f s , 9 — 9s+i, 9 — gs, then the relations (9.11), ( p. 12 ) form a special case of the 
difference Painleve IV equation (d— Pjv) of |L7]], §7. The parameters A, ao, ai, a 2 , 03 in our case are as follows: 

ao = — s — 2, ai = 1, a 2 = k, 

03 = s + 2 — k, A = ai + a 2 + a 3 + a = 1. 

Proof of Theorem \9. j . For the first part of the proof we do not need to distinguish between the cases where £ 7^ 
and £ = 0. Taking the residues of both sides of ([O]) at £ = s + 1 yields 

A s • ((s + 1)A + C,) = ((s + 1)A + C s+ i) • A s+1 , 



(fc + M 



— 1 —a s (3 s 
!/(<*.&) 1 



s + 1 + 6 S 6 S /3 S 
(r - £6 s )/& £(s + 1) + t - £6 S 



s + 1 + b 



s+l 



1& 



-1 -a s+ i/3 s+ i 
l/(a s+ i/3 s+ i) 1 



(fc + bs+l) A (r - a+i)//3 s+ i £(s + 1) + r - £6 S+1 
Comparing the diagonal terms on both sides, we get a system of two equations: 

(k + b s ) • [-(r - £o s )a s - (s + 1) - b s ] = (k + b s+1 ) ■ [b s+1 /a s+1 - (s + 1) - b s+1 



(9.14) 



(k + b s ) ■ [b s /a s + £(s + 1) + r - £6 S ] = (fc + 6 s +i) • [-(t - £b s+1 )a s+1 + £(s + 1) + r - £6 S+1 ] . (9.15) 

If we multiply (9.14) by a s +i, ( |9.15| ) by a s , and add the results, we obtain an equation which can be written as 
follows: 

(k + b s ) • [(1 - £a s )(l - a s+1 )b s + ra s (l - a s+1 ) + £(s + l)a s - (s + l)a s ] 

= (k + b s+ i) ■ [(1 - £a s )(l ~ a s +i)& s +i + ra 8 (l - a s+1 ) + £(s + l)a s - (s + l)a s ] . 

Now we subtract the LHS of the last equation from its RHS and divide the result by (1 — £a s )(l — a s+ i)(b s+ i — b s 
noting that {k + 6 s+ i)6 s+ i - (k + b s )b s — (k + b s+ i + b s )(b s+ i - b s ), we get 



(9.16) 



k + 0, 



s + l 



s + l 



s+l 



l-£a s 



1 - a 



= 0. 



s+l 



(9.17) 



Let us now ass ume that £ ^ 0. In this case, it is easy to se e that with the notation (9.6), the last equation is 
equivalent to (f/7|). To obtain (|J), we divide ( |9.15| ) by fl9.14|) , which yields 



b s /a s + t(s + l)+T -£b s 
-(T-£b s )a s -(s + l)-b s 



~(t - S,b s+1 )a s+1 + £(s + 1) + r - £6 s+ i 
fos+i/cts+i - (s + 1) - b s+ i 



(9.18) 



From ( J9.6p and (9.17), we have 

b s /a a + £(s+l)+T-£b s 



— (1 - £a s )6 s + £(s + 1) + r = L_^i( /s+1 - 1 - s ), 



-(t - £b s )a s -(s + l)-b s = - [(1 - £a s )fe s + ra s + s + l] = -(1 - £a s )/ s+ i, 
-(t - £6 s+ i)a s+ i + £(s + 1) + r - £6 S+1 = r(l - a s+1 ) -£•[(!- a s +i)o s+ i - (s + 1)] = (1 - a s+ i)(r + £/ s+1 + 
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■>s+l/ a s+l 



(s + 1) - b s+1 = — !— • [(1 - a s+1 )b s+1 - (s + l)a s+ i] = 1 : as+1 (s + 1 - / s+1 - k). 



a s +i a s +i 

This computation immediately implies tha t (|9.18|) is equivalent to ([T^) . To complete the proof of part (a) , we 
equate the (2,1) elements of both sides of ( |9.3|) , which gives 

(r - & s )/p s + (k + b s )/(a s f3 s ) = (r - + £(k + b s+1 )/(a s+1 s+1 ). 

We can rewrite the last equation as 

-j- [(r - £6 S ) + (A; + b,)/a.] =^—[(t- £b s+1 ) + £(fc + 6 s+ i)/a s+1 ] . 

Ps Ps+1 



It is easily seen to be equivalent to ( |9.9[ ), by (9.6). 
Now we assume that £ = 0. Then ( 9.17 ) becomes 



k + b 8 + b s+1 + ra s + s + 1 - 



s + 1 



= 0. 



(9.19) 



1 - a s+ i 

It is clear that with the notation ( |9.10 ), the last equation is equivalent to ( |9.11 ). To obtain ( 9.1 2| ) , we divide 
( |9.15D by ( |9.14| ) , which gives 

b s /a s +r 



-™ s +i 



From (|9.10|) and (|9.19|) , we have 



(s + 1) - 6 S 6 s+ i/a s+ i - (s + 1) - 



(9.20) 



'Js+i/a s +i 



(■s + 1) 



6 s /a s + r = — (tq s + 6 S ) = — (g s - s - 1), 

-ra s - 6 a - (s + 1) = -,g s , 
-ra s +i + t = (1 - a s+ i)r, 

[(1 - a s+ i)& s+ i - (s + l)a 8+ i] 



1 - a x 



a s +i 



-(s + l-k- g s ). 



This computatio n im mediately implies that (9.20) is equivalent to (9.12). Finally, we compare the (2, 1) elements 
of both sides of (3.3). This yields 



r k + b s 



(3. 



i.e., 



P. 



s+1 



s+1 



Ps 



a s 1 (ra s + k + b s )' 



which gives (9.13), completing the proof of part (b). 

10. A CONNECTION WITH THE q — P V I EQUATION OF M. JlMBO AND H. SAKAI 



□ 



10.1. Reduction to the q — Pvi system. In this section we show that the compatibility conditions for the Lax 
pairs corresponding to some of the families of polynomials orthogonal on g-lattices are equivalent to the q — Pvi 
system of M. Jimbo and H. Sakai (equations (19)-(20) in [l0|] ) for an appropriate choice of the parameters, or to 
a certain degeneration of this system. Thus, we now assume that the orthogonality set is of the form X = {q~ s }, 
where s runs either over Z> or over {0, . . . , N} (N 6 Z>o), and |g| ^ 0, 1. Hence we have a( — q( and r] = q. 
Then the corresponding Lax pair has the following form: 

X m s (C), m s (qO = M s (C)m s+ i(QD(C)-\ (10.1) 



I 



C-9- 



and the compatibility condition is (cf. equation (5.5)) 



qC - q 



M S (Q = M S+1 (C) [ I + 



(10.2) 
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To relate our situation to the one considered in pQ| , we make the following change of notation: x := Ci t :== Q 
Then we define 



A(x,t) = M s (x)((x-t)I + A s ), 
B {t) = -qtl-A^!, 
x(xI + B (t)) 



B(x,t) 



(x - qt) 2 



(10.3) 
(10.4) 

(10.5) 



Theorem 10.1. (a) The compatibility condition (10. i.) with s replaced by s — 1 is equivalent to the following 
equation: 



A(x,qt)B(x,t) = B{qx,t)A{x,t). 



(10.6) 



(b) If the matrix D(Q — diag(di(£), <fe(0) * s linear in (, so that di(Q = X\(C — a^), (^(C) = Aa(C — a i) with 
Ai , A2 7^ 0, then the matrix A(x,t) is quadratic in x, and if we write 

A(x, t) = A (t) + Ax{t)x + A 2 x 2 , 



have: 



A 



I , ki = q Ai, k 2 = g fc A 2 

Ao(f) has eigenvalues t#i, £6*2, 
det A(x, t) — kiK2(x — t) 2 (x — 0,3) (x — a&), 



(10.7) 

(10.8) 

(10.9) 
(10.10) 



and the parameters kj, dj, 9j are independent oft. 



Recall that the natural number k in (10.8) defines the asymptotics of the solutions m s (C) as £ — > 00 (see 
Theorem O). 



Remark 10.2. The equation ( 10. q ) can be viewed as the compatibility condition for the following pair of q- 
difference matrix equations, cf |L0|]: 

n(x, qt) = B (x , t)n(x , t) , n(qx,t) = A(x,t)n{x,t). (10.11) 

One way to construct a matrix n{x, t) which solves the system ( |l0.1l| ) is as follows. Let V(() = diag(ui(£), ^(0) 
be a diagonal matrix such that V(qQ = D(()V(() for all £. Then for all s, we define a new matrix n s (£) by 

' c s q {l) nt-oofe'c - V' 1 • • ncr 1 if m > i, 

C s <7® nt~(9*C - 1) • «.(C) • ^(C)- 1 if < M < 1. 



m s (C) 



(10.12) 



Substituting this into the first equation of the Lax pair (10.1), replacing s by s — 1 and simplifying yields 



C • n s (C) = ((C - <T S+1 )/ + • n.-i(C). 



Since j4f_j = 0, this is equivalent to 



C(q + (-g- 8+1 J-^-i)) ... 

^-i(0 = — 77 rTTTTo -risiO- 



(10.13) 



(C-Q- s+1 ) 2 

On the other hand, if we substitute the first equation of the Lax pair ( lO.lj ) into the second one, use (10.12) and 
simplify the result, we obtain 



n s (qQ = M S (C) ■ ((C - Q- s ) + A.)n.(Q. 



(10.14) 



Now if we l et x = (, t = q ~ s and n( x t) = n s ((), then with the notation ( |10.3[ ), ( |10.4| ), ( |10.5| ) of Theorem |l0T 
the system (110.131), (110.141) leads to fllOTlJ). 
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Proof of Theorem 10.1. (a) Since A 2 S _ 1 = 0, we have 



A(x, qt)B(x, t) = M a _ x (C) • ((C - q- s+1 )I + A s _!) 



C((C - q- s+1 )i - As-!) 



(C 



CM s -i(C) 



and 



B(qx,t)A(x,t) 



q((q(( - q- s )I - A s -!) 



ga(C-5-) a 

s— 1 



C 



M s (0-((C-g- s )i + A s ) 
A, 



M S (C)- 1 + 



C-9- 



Hence, if we multiply both sides of ( 10. 6| ) by (j — (£ — q 
by £, we obtain 



C-q- 



l q- X A s -i 
M a _i(C) = M a (C)- 



r J = ^+(c- 



g .A a _i and divide 



-4., 



Replacing s by s + 1 yields (10.2). 

(b) Since the matrix D(C) = diag(Ai(£ — 03), A2(C — 0,4)) is linear in £, it follows from Theorem |IP|(c) that we 
can write 

K 2 ( 



Ms(C) 



C s 



for a constant matrix C s , where ki,k 2 are given by (10.8). This immediately implies that the matrix A(x,t) is 



quadratic in x with the leading coefheient A2 = diag(Ki, K2). Since det M S {C) = det D(Q for all £ by Lemma 2.1 
and the matrix A s is nilpotent, we see from (10.3) that det A(x, t) — kik 2 (:z; — t) 2 (x — a^)(x — 04). In particular, 
taking x = 0, we see that detAo(t) = t 2 5, where 6 is a constant independent of t. To complete the proof, it 
therefore suffices to show that Tr Ao(t) = where r is a constant independent of t. To this end, we study the 
compatibility condition ( 10. 2| ). Taking residues of both sides of (10.2) at £ = q^ 11 ^ 1 yields 



As- 



Kiq 





--S-1 



Cs = q ■ 



Kiq 







K 2 q 

If we compare the asymptotics of both sides of ( 10.2] ) as £ — > 00, we get 





K 2 g 





-s-l 



Cs 



A, 



C s + A s 



Kiq 




Kl 
K 2 



.4 



s+l- 



Multiplying (10.16) by q s and subtracting it from (10.15) gives 

(A. - q- s I)C s = qC s+1 (A s+1 - q- s - l I). 
Since A (t) = C S (A S - q~ s I), this shows that Tr A (qt) = q • Tr A (t) for all t. 



(10.15) 



(10.16) 



□ 



10.2. Solution of the q — Pvi system. We are now in a position to quote a result of M. Jimbo and H. Sakai 
[ p7o| . The situation considered in their work is more general: it is assumed that B(x,t) has the form 

x(xI + B (t)) 



B(x,t) 



(x — qta\){x — qta,2) 



and instead of ( 10.10 ), it is assumed that 

det^4(x,i) = kik 2 (x — ta\)(x — ta,2)(x — a^)(x — 04). 



(10.17) 
(10.18) 



Thus, the situation of Theorem 10.1 corresponds to the specialization a\ = a 2 = 1. Now we have 



Theorem 10.3 (Jimbo-Sakai, fllO[). Assume (10.17), (10.",), (10.6), (10. i) and (10. IS). Also, suppose that 



yj, w } 



^0 (7 = 1,2), Oj ^0 (j = 1,2,3,4) and k 1 ^k 2 . 



:!0 
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(a) Define y — y(t) and Zi = Zi(t) (i = 1,2) by 

Ai2(y,t) = 0, A n (y,t) = K1Z1, A 22 (y,t) = k 2 z 2 , 

where are the elements of the matrix A(x,t), so that z\z 2 = (y — ta{)(y — ta 2 )(y — a 3 )(y — a 4 ). In 

terms of y, Z\, z%, the matrix A(x,t) can be parameterized as follows: 



where 



a 



A(x,t) 
1 

Kl - K 2 
1 



Kl - K 2 



_ / ki ((x - y)(x - a) + zi) k 2 w(x - y) 

\ kiw^ 1 (jx + S) k 2 ((x - y)(x - (3) + z 2 ) 

V 1 ((^1 +02)t — K\Z\ - k 2 z 2 ) - k 2 ((ai + a 2 )t + a 3 + a 4 - 2y) 
-y^ 1 ((#1 + 9 2 )t - K1Z1 ~ k 2 z 2 ) + ki ((ai + a 2 )t + a 3 + a 4 - 2y) 



7 = zi + z 2 + (y + a)(y + f3) + (a + (3)y - aia 2 t 2 - (ai + a 2 )(a 3 + a 4 )i - a 3 a 4 , 
(5 = y^ 1 (aia 2 a^a A t 2 - (ay + zi){j3y + z 2 )). 



(b) Define z = z(t) by 



(y - tai)(y - ta 2 ) 



qKiz 



z 2 = qni{y - a 3 )(y - a^z. 



Introduce the notation y — y(qt), z — z(qt), w = w(qt), and set 



a x a 2 a x a 2 
°i = — n — , o 2 = — — , 6 3 



1 

qui 



1 

K 2 ' 



Then the compatibility condition (10. t) is equivalent to the following system of equations: 

yy _ (g ~ tbi)(z - tb 2 ) 
a 3 a 4 (z - b 3 )(z - b 4 ) 
zz _ (y - tai)(y - ta 2 ) 
6 3 6 4 (y - a 3 )(y - a 4 ) 
w 6 4 z — 63 
w 6 3 z — 6 4 



(10.19) 
(10.20) 
(10.21) 



Remark 10.4. Equations ( 10.1 9| ) — ( |l . 2 1 ) allow us to compute (y,z,w) if we know (y,z,w), and vice versa. 



Remark 10.5. Unfortunately, it was beyond our technical abilities to follow the proof of Theorem 10.3 in |10| . 
However, we were able to verify the statement of the theorem using computer simulations with random values of 
the parameters Kj , aj , 9j . 



10.3. Degeneration of q — Pvi- We have mentioned in subsection 7.5 that the compatibility conditions for Lax 
pairs corresponding to certain families of orthogonal polynomials are equivalent not to the q — Pyi system but 
to a degeneration of it. We now describe this degeneration. 



Theorem 10.6. With the notation of Theorem 10.1, suppose that the matrix D(C,) = diag(di(£), cfe(C)) * s su ch 
that di(C) = Ai(C — a 3 ), cfe(C) — ^2, where Xj,a 3 ^ 0. Then 
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(a) The matrix A(x, t) is quadratic in x, and if we write 

A{x, t) = A (t) + Ax{t)x + A 2 x 2 , 

we have: 

K i 1 K° — n k \, 
/ ' 1 — 



A- ? 



A (t) has eigenvalues t0\, t9 2 , 
det A(x, t) = k°k 2 (x — t) 2 (x — a 3 ), k 2 



<r fe A 2 , 



and i/ie parameters k°, 6°, a 3 are independent oft. 
(b) We can parameterize the matrix A{x,t) as follows: 



a(t _ ( Kl({x-y°){x-a°) + zl) n° 2 w°(x-y°) \ 



where 



1 r 



(y )- 1 ^! + d° 2 )t - k\zI - K ° 2 z°) + K ° 



(10.22) 

(10.23) 

(10.24) 
(10.25) 

(10.26) 



O O / O i O \ O i C\j_ i o 

7 =^2-(y +a)-y +2t + a 3 , 
S° = (y°y 1 (-a° 3 t 2 + (a°y° + z°)(y° - z 2 )). 



Define z° = z°(t) by 



(y°-t) 2 

qK°z° 



O O / O \ C 

z 2 = QKiiV - a 3) z 



Introduce the notation y° — y°(qt), z° — z°(qt), w° = w°(qt), and set 

7 O _ 7 O _ 7 



Then the compatibility condition (10. t) is equivalent to the following system of equations: 

y°y° (z° -ibl){z° -tb° 2 ) 



t, 2 u, 3 



(y°-t) 2 



w 
w 



(10.27) 
(10.28) 
(10.29) 



Proof, (a) The proof of this part is almost identical to that of Theorem 10.1(b) and will therefore be omitted, 
(b) One checks directly that all formulas of [fol are compatible with the following limit transition: 



Ki — » /%, ai, a 2 



6 2 ^6. 



2 , a 3 ->a 3 , a 



y -> V , Z\ 



K 2 — > 0, K 2 Q4 



K 2 /3 



o o c o ro 

K 2 J — > K 2 7 , K 2 d — > K 2 d , K 2 ^2 



o o o o 

t 2 Z 2 , K2W — > K 2 W . 



This li mit transition takes the parame te rizatio n of A(x,t) gi ven in T heorem 10.3 (a) to the parameterization 
( |10.26|) , and it takes the system (|l0 19|) - (|l0.2l|) to the system ( |l0.27|) -( |ToT2g| ) . □ 
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11. Applications: recurrence relations for some polynomials of the Askey scheme 



11.1. Notation. In this section we illustrate the results of §§| - 10 by considering several specific examples: Char- 
lier polynomials, Meixner polynomials, Krawtchouk polynomials, g-Charlier polynomials, little g-Laguerre/Wall 
polynomials, alternative q-Charlier polynomials, little g-Jacobi polynomials and q-Krawtchouk polynomials. In 
the first four cases we solve the compatibility condition explicitly, we write down a recurrence relation for the 
corresponding Fredholm determinants in terms of the solution, and we provide the initial conditions for all of 
our recurrence relations. This is done in subsections [11. 2\ |1 1 .4 , 11.5 and 11.6. For the other four families, we 
contend ourselves with making some general remarks in subsection 11.7. 

The general notation of this section is that of §[| and §|| Recall that we are considering a family {P n (0}n=o 
of monic polynomials orthogonal on a discrete, but not necessarily locally finite subset X = {tt x } x=0 of K (where 
N e Z>o U {oo}), with respect to a strictly positive weight function uj : X — > K. If k is a natural number, k < N, 
we can consider a kernel K on X x X defined by the formula 



\\P k ^\\Z 2 uj(x)(Pk^x)Pk-iM - PLxi^PM), x = y, 



(11.1) 



1/2 

where ||-Pfc-i|| w = [Pk-i,Pk-x)J denotes the norm of Pk-i(() with respect to the inner product defined by u>. 
Up to conjugation, this coincides with the kernel introduced in the beginning of §|| (see equation (4.1)). For all 
k < s < N, we define a subset 2) s = {tt x } x=s Q X, and we are interested in the Fredholm determinants 



D, = detfl - K 



(11.2) 



11.2. Charlier polynomials (jT3|, §1.12). The n-th Charlier polynomial is defined by 

C n (x;a) = 2 Fo[ 
These polynomials satisfy the orthogonality relation 

oo x 

^2 ~ C m (x; a)C n (x; a) = a~ n e a n\S mn , 



(11.3) 



x=0 



(11.4) 



where a > 0. Thus the orthogonality set for Charlier polynomials is X = Z>o, and the weight function is 
lu(x) = ^j. The polynomial C n {x\ a) is not monic in general; in fact, its leading coefficient is (— a)~ n . Hence the 
corresponding family of orthogonal polynomials (recall that the orthogonal polynomials that we use are monic, 
see §||) is {P n (C) — (~ a ) n Cn(C; a)}^ =0 . We call P n (C) the n-th normalized Charlier polynomial. Now from (11.4), 
we find that (P n , P n )u> = a n e a n\ for all n > 0. After these preliminaries, we can state our main result for Charlier 
polynomials. 



Theorem 11.1. If K is the kernel corresponding to the family {-Pn(C)}^o of normalized Charlier poly- 

nomials, then the Fredholm determinants D s defined by (11. i) can be computed from the following recurrence 
relation: 



s+2 
s+1 



D 



s+1 



f; 



(s+iy. 



(9s 



(11.5) 
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Here, the scalar sequences {e s } s >k, {fs}s>k, {gs}s>k and {h s } s >k satisfy the following recurrence relations 

ae s 



e s +i 

fs+l = 



f s (g s + k-s-l) 
ag s 



f s (g s 



9s+i 



-«-l)Cff. 
s + 1 



+ k-s-iy 

g s — k + 2s + 3, 
h s . 



fs+l 1 - fs+l 
h s+ i = a^ 1 ■ f s ■ (g s - s - 1 

The initial conditions for the recurrence relations fill.!\ )- ftll.£> ) are given by 

D k+1 =e- ak ■^(-k-l--a), 

a k ■ (fc-1)! 



D, 



e-k 



fk — - 



$(l-fc;l;-a)' 
a-$(l -k;2;-a) 
$(l-jfc;l;-a) : 



9k 



(k + 1) ■ $(1 - fc; 1; -a) • $(-fc; 2; -a) 
$(-fc;l;-a)-$(l-fc;2;-a) : 



(11.6) 
(11.7) 

(11.8) 
(11.9) 

(11.10) 
(11.11) 

(11.12) 

(11.13) 
(11.14) 



Remark 11.2. As we have have seen in 
special case of the d — Pjy equation of \ 



$(u;w;z) = iFi\ 

9 (cf. Remark 9.4(b)), the recurrence relations < $Tj\) , JlL^) form a 



Proof. The proof is a straightforward application of our previous results. For the reader's convenience we provide 
some remarks; the same ones apply to Theorems 11.4 and 11 .6| , and hence the proofs of those results will be 
omitted. 

To make the notation of the present section more uniform, we have been writing e s for [3 S and h s for m s , 
where {(3 S } and {"i^ 1 } are the scalar sequences defined in §|| andj P , resp ectively. The symbols f s and g s have 
the same meaning a s in §|^. Then the recurrence relations ( |ll.E| )- (|ll.9| ) are obtained directly from Theorem 
l^b) and Theorem j.3(b) . T o find the initial conditions, one uses the definitions of bk, ctk, Pk, fk, 9k, tog ethe r 
with Propositions 3.1, 3.2, 6.6, and the obvious identities h^ = rnj 
&■ 



k\, bkPk = —qk, which follow from (|6.1 



□ 



11.3. We now illustrate the concluding remark of §|3J by showing how one can use Proposition 3.3 to obtain 
difference equations satisfied by orthogonal polynomials. 

Proposition 11.3 (cf. jl^], §1.12, equation (1.12.5)). The k-th normalized Charlier polynomial Pk(C) solves the 
following difference equation: 

-kP k (0 = aP k (( + 1) - (C + a)P k (C) + (Pk(C - I)- (11-15) 

Proof. We use the notation of Proposition [3.3[ Recall that in the case of Charlier polynomials, we have D(() = 
diag(C, a). We also observe that from the proof of Theorem [2.4| it follows that the matrix mx{C) has a full 
asymptotic expansion in £ as C —> oo; in particular, we can write, by (|2.4|), 



mx(C) 



C k o 
o c k 



a (3 
7 5 
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Therefore, as £ 
M(C) 



oo, we have 

= m x ^-l)-D(Q-m^(0 



m x (C-l)- 
m x (C) ■ 



o 



(c-i; 



C~ k o 
o c fc 



1 + aC" 1 /?C _1 

7 C _1 i + ^C 1 

7 a 



C o 

a 


(C - i) fe 

-1 



C-fc o 

a 



C fc 

o c k 



■m x \C) 



(i-i/C) k o 

o (i - i/cr k 



C o 

a 



1 - aC 1 -PC 1 
-iC 1 i-SC 1 



+ 0(C 1 ) 



)+o(r 1 ). 



Since M(£) is entire by Prop o sitio n 3.3, the last term 0(£ 1 ) is identically zero by Liouvillc's theorem. Hence 
the system of equations (3J3), (3^) takes the following form: 

C-P fc (C-i) = (C-*0-P fe (C)-/3-cP fc -i(O, (n.16) 

C-cP fe _!(C-l) = 7 -Pfc(C)+a-cP fe -i(C). (11.17) 

Note that since det M(£) = det -D(£) for all ( by Lemma |2.1| , we have f3j = ak; in particular, (3 ^ 0. Now from 
( |11.16[ ), we find that 

cPfe-i(C) = - & ■ [(C - *)flfe(0 - CPfe(C - 1)] ■ 
Substituting this into ( |11.17| ), multiplying the result by (3 and using /37 = ak, we obtain 

C • [(C - 1 - fc)P fc (C - 1) - (C - l)Pfc(C - 2)] = ak ■ P k (0 + a ■ [(C - fc)P fe (C) - CPfc(C - 1)] • 
Dividing the last equation by ( and replacing £ by C + 1> we find that it is equivalent to □ 

11.4. Meixner polynomials §1-9). The n-t/i Meixner polynomial is defined by 

1' 



M n (x;/3,c) = aFi 
These polynomials satisfy the orthogonality relation 



—n, —a; 




1 - 



E 



^-c x M m {x;(3,c)M n {x;(i,c) 



(11.18) 



(11.19) 



where /3 > and < c < 1. Thus the orthogonality set for Meixner polynomials is X = Z>o, and the weight 
function is uj{x) = ^rc x . The leading coefficient of the polynomial M n (x; (i, c) is ^j^y^—, so the corresponding 
family of monic orthogonal polynomials is {P n (C) = (/3) n (l — l/ c )~™-<Wn(C; &> c )}^Lo- We call P n (C) the n-i/i 



normalized Meixner polynomial. Now from ( 11.19 ), we find that (P„,P„) 
have 



(l-c) 



for all n > 0. Then we 



Theorem 11.4. //if is the kernel corresponding to the family {Pn(C)}5^Lo of normalized Meixner poly- 

nomials, then the Fredholm determinants D s defined by (11. i) can be computed from the following recurrence 
relation: 



s+2 

S+l 



D 



s+l 



p + s ' (s + l)! 



1 + cg s 



[(1 + cg s )f s+ i 



1] -hi 



(11.20) 
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Here, the scalar sequences {e s } s >k, {fs}s>k, {gs}s>k and {h s } s >k satisfy the following recurrence relations 

ce s g s (1 + g s +i)f s +i + {P + k - l)g s +i -s-1 



e s +i 



9s+i 

f s +i = l - P - k - f t 



(1 + cg s )f s+1 +k-s-l 
s (3 + s 



9s+i 



1 + g s 1 + cg s 
(/ s+1 -l-s)(/ s+1 -l-s + fc) 

cg s fs+i{fs+i + P + k - l) 



_ (1 + cg s )(s + 1 - f s+1 ) 

n-s+l — 7-5— — r • Us- 

c(p + s)g s 

The initial conditions for the recurrence relations fill.2(\) -( 11.2$) are given by 

D k = (1 - C )^-D, Dk+1 = . c k . (1 _ c) K0+k-x) . F( _ fc; _ k p 1/c) 



efc 



(3c- (k- 1) 



|2 



fic 



F(l-fc,l-fc;l + /3;l/c)' 
/fc = 0, 

fc F(l-fc,l-fc;l + /3;l/c) 



F(-fc,l-fc;/3;l/c) 



(11.21) 
(11.22) 
(11.23) 
(11.24) 

(11.25) 

(11.26) 
(11.27) 
(11.28) 
(11.29) 



Proof. See the proof of Theorem 



F(u,v;w;z) = 2 Fi 



u, v 
w 



11.1 



□ 



Remark 11.5. As we have have seen in §0 (cf. Remark |9.4| (a)), the recurrence relations ( 11.22| ), ( 11.23| ) form a 
special case of the d — Py equation of [|l7|. 



11.5. Krawtchouk polynomials ( |[L3| , §1.10). The n-th Krawtchouk polynomial is defined by 

r 



K n {x;p,N) = 2F1 
These polynomials satisfy the orthogonality relation 



-n, —x 
-N 



P 



N 

E 

x=0 



N 



p x (l - p) N - x K m (x;p, N)K n (x; p, N) 



(-l)"n! (I -p\ n 

(-N)n 



(11.30) 



(11.31) 



where N 6 Z>o and < p < 1. Thus the orthogonality set for Krawtchouk polynomials is X = {0, ... , N}, 
and the weight function is ui{x) — {T)p x {i — p) N ~ x . The leading coefficient of the polynomial K n (x;p, N) is 
(— TV),^ 1 ^", so the corresponding family of monic orthogonal polynomials is {P n (Q — (—N) n p n K n ((;p, N)}^ =0 . 
We call -Pn(C) the n-th normalized Krawtchouk polynomial. Now from (11.31), we find that {P n ,Pn)uj = 
(-l) n nl(-N) n p n (l - p) n for all < n < N. Then we have 

Theorem 11.6. If K is the kernel corresponding t o the family {P n (C)}n=o °f normalized Krawtchouk 

polynomials, then the Fredholm determinants D s defined by (ll.i) can be computed from the following recurrence 
relation: 

p s-l(l-p)N-s+l l+pg B /(p-l) 



D, 

TJ. 



+2 
+1 



D 



s+l 



N 
s + l 



(N - s) 2 



e s gi 



(l+pg s /(p-l))f s+1 -s-l 



(11.32) 
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Here, the scalar sequences {e s } s >k, {fs}s>k, {9s}s>k and {h s } s >k satisfy the following recurrence relations 

pe s g s (1 + sf 8 +i)/s+i + (k-N - l)g B +i -8-1 



e s +i = 



(1 -P)9s+i (l+pg s /(p- l))/ s +i +k-s-l 

s (l-p)(N-s) 



f s+ i = N + l — k — f l 
9s+i 



l + g s p-l+pg s 
(l-p)(f s+1 -l-s)(f s+1 -l-s + k) 



P9sfs+i(N + l-k- fs+i) 
, (p - 1 + pg s )(f s+1 - s - 1) 

«s+l — 7TT ^ ' Us 



p(N - s)g s 

The initial conditions for the recurrence relations (11. Si )-( !!. 3oQ are given by 

D k = (l- p)*C"+i-*) , D k+1 = ■ p k ■ (1 - p)^" fc ) • F(-k, -k; —N; 1 - 1/p), 



Npil-p^^ik-iy. 2 



9k 



F(l-k,l-k;l-N;l-l/p) 

h = 0, 

k(l-p) F(l-k,l-k;l-N;l-l/p) 
Np F(-k,l-k;-N;l-l/p) 
h k = k\, 



(11.33) 
(11.34) 
(11.35) 
(11.36) 

(11.37) 

(11.38) 
(11.39) 
(11.40) 
(11.41) 



where 



F{u, v, w; z) = 2-F1 



u, v 
w 



Proof. See the proof of Theorem 11.1 



□ 



Remark 11.7. As we have have seen in Jp (cf. Remark |9.4| (a)), the recurrence relations (11.34), (11.35) form a 
special case of the d — Py equation of JL7|. 



11.6. g-Charlier polynomials (|13|, §3.23). In this subsection, we assume that q is a fixed real number, 
< q < l.Thc n-th q-Charlier polynomial is defined by 

, ,, .r»n+l' 

C„(C;a;g) = 2^1 







q" 



These polynomials satisfy the orthogonality relation 



Y] , a x q^C m (q x ;a;q)C n (q x ;a;q) = q n (-a;q) 00 (-a 1 q 7 q;q) n 5 ri 



x=0 



(11.42) 



(11.43) 



where a > and 



(-«; 9)00 = + a <?')- 



1=0 

Thus the orthogonality set for q-Charlier polynomials is X = {q~ x }^L , an( i the weight function is lo(x) — 

{q-q) Q^ 2 ^ ■ T ne lading coefficient of the polynomial C„(£;a;q) is (— i) n q n a~ n , so the corresponding family of 

orthogonal polynomials is \P n (Q = (—l) n a n q~ n C n (C,; a; <?)}^L - We call -Pn(C) the n-th normalized q-Charlier 
polynomial. Now from ( 11.43 ), we find that 

(P n , Pnh = a 2n q- 2n2 - n (-a; qUi-a^q, q; q) n 
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for all n > 0. Then we have 

Theorem 11.8. If K is the kernel \ll.\ ) corresponding to the family {-Pn(C)}5^=o °f normalized q-Charlier 
polynomials, then the Fredholm determinants D s defined by (11. i) satisfy the following recurrence relation: 

(11.44) 



Ds+2 

D s +i 



D 



s+l 



where 



■qC* 1 ) -u s -h 2 s , 



u s = q " ■ ■ {PsPs + aq q s ) 

q s 



(11.45) 



for all s > 0. The scalar sequences {p s } s >k, {qs}s>k, {Ps}s>k cind {/i s }s>fc can be computed from the following 
recurrence relations (which involve additional sequences): 

e s = a(q- s - 1 - 1) + q k - 1 (aq~ k p s - r s (3 s ), 

(p s (3 s + aq~ k q s ) ■ (r s a s - p s j s + g fc-a_1 r a ), 

q s +i = q^q^e' 1 ■ (p s (3 s + aq~ k q s ) 2 , 



Ps+1 



q ; /<. V ; 



r s+ x =q r s e 



■'- 1 1 (r.a.-jvy. + g*— V,) 2 , 



a s+1 = a s + q k ~ 1 p s - q k p s +l, 
[3 S+1 =p„- q k q s+ i, 

7s+i = 7s + q k ~~ 1 r s , 

_ p s f3 s + aq~ k q s 

"s+i — • n>s' 

aq s 



The initial conditions for the recurrence relations ( 11.44 ) an d ( H-V' )-( H.5$() are provided by 

fe-i 

D k = (-a; q)-J ■ o-(5) ■ [] [(-a- X q; qVq^} , 



ra=0 



A— 1 



D 



fe+i 



(q;q)k 



■ q~® ■ G q (q- k , q- k ; -q 2k /a) ■ [J [(-a^q; g^gO^ 



71=0 



Pk = (1 - 9 



G g (g fc ,g k ;-q 2k /a) 
qu = (q; q)l • <T fc(fc+1) • G,(g- fe , <T fc ; -^/a)" 1 



g fc 2 (l_ g -fc) G g (g- fe , g 1 - fe ;-(? 2 ' £ - 1 /a) 2 



(q;q)k(q;q)k-i G q (q 



fe /c ■ 2 & 



q 2k /a) 



where 



a k = -l-q p k , 
Pk = -q k qk, 

k 2 -l 

1 n f„l-k „\~k. „2k-2 /„\ 

7*; = 7 72 G q (q ,q ;-q /a), 

frfc = <T fc2 • 

G q (u,v;z) = ■ 2 fj) [ : 



(11.46) 
(11.47) 
(11.48) 
(11.49) 
(11.50) 
(11.51) 
(11.52) 

(11.53) 



(11.54) 

(11.55) 

(11.56) 

(11.57) 

(11.58) 

(11.59) 
(11.60) 

(11.61) 

(11.62) 
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Proof. As before, the proof is quite straightforward. We have been writing h s for m* 1 (defined in §||). The 
formulas (11.45) and ( 11.46|)- ( 11.52D follow immediately from Theorem S.2 (note that S s = aq~ k for all s, since 
we have 5 s +i = S s from (8.11), and 6k = aq~ k from ( 5.1 C|) ) . Then (11.44) and (11.53) are deduced from (4.7) 
and ( |4.6|), r espectively. Finally, the initial conditions ( 11.54] )— (11.62) are easily obtained from Proposition 6.6, 
"T2| and 6.51. □ 



Remark 11.9. As we have mentioned in subsection 7.5, the recurrence relation (11.47)— ( 11.52) for g-Charlier 
polynomials is in fact equivalent to a certain degeneration of the q — Pyi equation of 1 1 . This is a special case 
of Theorem 10.6. 



Remark 11.10. In the case of g-Charlier polynomials, it is possible to solve the compatibility condition for the 
corresponding Lax pair by a method similar to the one used in Theorem |9.3| (b) . Namely, it is easy to see that 
the Lax pair can be parameterized as follows: 



m s+ i 



K' = { /+ «-'"'>"'(-p1, 



p s d s C s 
Cs) -Ps 



m (at)- ( 9 k (C-l) + b s b s c s \ / ( C - l)" 1 

m *WO-y aq -k /Cs aq - k ) ms+l[Q { a- 1 

Then the compatibility condition gives the following recurrence relations for the parameters a Sl b s , c s ,p s 

p s (b s + aq- k a s )(q k - s - 1 - q k + b s + aq- k a s ) 



Ps+1 = 



q k p s (b s + aq k a s ) + q ■ (q k s 1 - q k ) ■ aq k a s 



b s+ i =b s + q k x p s - q k p s +i, 



q^Ps+i + aq n a s 



- q k + b s - q k p s +i + aq k a s 



c s +i 



aq k a s c s 



aq k a s - q k 1 p s 



With this notation, the recurrence relation for the Fredholm determinants is ( 11.44 ), the same as before, but 
now we have 



q k - (aq k + b s /a s ) 
a s c s 



and the recurrence relation for h s is given by 

h s+ i = a^ 1 ■ (aq~ k + b s /a s ) ■ h s 



The initial values ak,bk,Ck can be easily found from ( 11.56 ) — (p.1.61 ). 

The main difference between this situation and that of §|9| is that we cannot further reduce the recurrence 
relations for a s ,b s ,p s to relations involving only two sequences of parameters. So this method cannot be used 
to show that our recurrence relations in the case of g-Charlier polynomials are equivalent to one of H. Sakai's 
g-difference equations |]l7| . From the computational point of view, this method is slightly easier to use than the 
one presented in Theorem U.S. 
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11.7. Concluding remarks. It follows from Theorem 10.1 that the recurrence relations corresponding to the 
little g-Jacobi polynomials and the g-Krawtchouk polynomials can be reduced to special cases of the q — Pyi 
system of |hJ. Also, it follows from Theorem 10.6 that the recurrence relations corresponding to the g-Charlier 
polynomials and the little g-Laguerre/Wall polynomials can be reduced to special cases of a certain degeneration 



of the q — Pvi system described in subsection 10.3. However, it is more convenient to use the formulas of §|8| for 
practical computations. In addition, the method of §|| covers the case of the alternative g-Charlier polynomials, 
whereas we do not know if the recurrence relation corresponding to these polynomials can be reduced to one of 
the equations of H. Sakai's hierarchy. 

As far as using the formulas of §|| is concerned, there is no essential difference between the g-Charlier poly- 
nomials and the other four families of basic hypergeometric orthogonal polynomials that we consider here. So 
we have decided not to write out explicitly the results we have obtained for these four families. On the other 
hand, we have carried out all the calculations for some specific values of parameters in Maple, and in § |l2] we 
present a few plots of the "density function" (difference or g-derivative of D s ) for the eight families of orthogonal 
polynomials considered in this section. 



12. Numerical computations 



12.1. The plots in this section have been obtained in Maple by using the formulas of 9q and subsections 11.2 



11.4 and 11.5 for the specific values of parameters indicated below. 



12.2. The following are two plots of the density function D s+ i — D s for the family of Meixner polynomials. The 
parameters (cf. subsection 11.4 or subsection 7.2) are k — 4, c = 0.01, (3 = 3000 for the first graph and k = 4, 
c = 0.9, (3 = 0.5 for the second graph. The ^-coordinate in each case is s. 



12.3. The following are the plots of the density function for the famili es of Charlier, g-Charlier and alternative 
g-Charlier polynomials (left to right). The parameters (cf. subsection 7.2) are k — 6, a — 20 for the Charlier 
polynomials (first graph) and k = 6, a — 20, q — 0.96 for the g-Charlier and the alternative g-Charlier polynomials 
(last two graphs). In the case of Charlier polynomials we plot the difference derivative, D 8 +\ — D s , of D s , and 
the x-coordinate is s, while in the other two cases we plot the g-derivative, q s ■ (D s+ i — D s )/(1 — q), of D s , and 
the ^-coordinate is q~ s . 
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12.4. The following are the plots of the density function q~ s ■ (-D s +i — D s )/{\ — q) for the families of little 
g-Laguerre /Wa ll polynomials (first graph) and little g-Jacobi polynomials (second graph). The parameters (cf. 
subsection 7.2) are k — 6, a = 0.5, q — 0.9 for the little q-Laguerre polynomials and k = 6, a = 0.5, b = 1.5, 
q = 0.9 for the little g-Jacobi polynomials. The ^-coordinate in each case is q s . 



12.5. The following are the plots of the density function for the families of Krawt chou k polynomials (first 
graph) and g-Kratwchouk polynomials (second graph). The parameters (cf. subsection 7J2) are k = 5, N = 80, 
p = 1/(0.7 + 1) for the Krawtchouk polynomials and k — 5, N — 80, p = 0.7, q — 0.98 for the q-Krawtchouk 
polynomials. In the first case we plot the difference derivative, Z? s +i — D Sl of D s , and the x-coordinate is s, while 
in the second case we plot the normalized g-derivative, (q~ N — 1) • q s ■ (D s+ i — D s )/(1 — q)/N, of D s , and the 
x-coordinate is N ■ (q~ s — l)/(q~ N — 1). 
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